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SIGNAL-TO-NOISE  RATIO  REQUIREMENTS  FOR  HALF-WAVE  AND 
FULL-WAVE  NONLINEAR  DETECTORS  WITH  ARBITRARY  POWER  LAWS, 
SAMPLING  RATES,  INPUT  SPECTRA,  AND  FILTER  CHARACTERISTICS 


INTRODUCTION 


The  purpose  of  this  report  is  to  determine  the  signai-to-noise  ratio 
requirements  for  various  half-wave  and  full-wave  rectifiers,  arbitrary  input 
spectra,  and  post-detector  filter  characteristics.  Both  continuous  and 
sampled  systems  are  considered,  as  well  as  broadband  and  narrowband  spectra. 

The  system  of  interest  is  indicated  in  figure  1.  The  input  x(t)  is  a 

2 

real  stationary  zero-mean  Gaussian  process  with  variance  a  and  normalized 
correlation  p(T) .  The  following  analysis  will  utilize  these  general 
parameters  where  possible;  however,  since  the  system  of  figure  1  is  to  be  used 
for  detection  or  a  decision  between  two  hypotheses,  we  will  later  specialize 
to  the  cases 


(1) 


where  the  subscript  denotes  the  hypothesis  number,  \  or  0.  That  is,  H 


denotes  the  signal  plus  noise  hypothesis,  while  Hq  denotes 
hypothesis. 


the  noise-only 


1 


X  (  t  ) 

Non  1  inear  f  ty 

y(t) 

Fi Iter 

z  ( t ) 

gCx) 

h(T),HCf ) 

Figure  1.  Detection  System 


Figure  2.  Discrete  Detection  System 


—  w  — 

ffj  —  w— * 

2W 

1_ i 

f 

-fc 

0 

Figure  3.  Flat  Bandpass  Input  Spectrum 
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The  nonlinearity  g{xl  is  characterized  by  output 


y(t)  =  gfx(t)]  ; 


(2) 


this  is  a  memoryless  transformation  of  the  input  at  the  same  time  instant.  We 
shall  be  particularly  interested  in  the  class  of  half-wave  rectifiers 


(3) 


and  the  class  of  symmetric  full-wave  rectifiers 


g{x]  =  1x1  ”  for  al  1  X  . 


(4) 


The  particular  rectifier  with  «  =  0  is  not  useful  for  detection  purposes  here, 
in  either  the  full-wave  or  half-wave  case.  For  the  full-wave  rectifier,  u  =  0 
corresponds  to  a  constant  output,  y(t)  =  1,  regardless  of  what  the  input  x(t) 
is.  For  the  half-wave  rectifier,  w  =  0  corresponds  to  y(t)  =  1  whenever 
x(t)  >  0;  but  if  we  are  trying  to  decide  between  two  zero-mean  processes  of 
different  levels,  this  information  is  lost  at  the  half-wave  rectifier  output. 
Hence,  we  assume  w  >  0  from  this  point  on,  when  we  deal  with  u-th  law 
rectifiers . 

The  low-pass  filter  is  characterized  either  by  its  impulse  response  h(T) 
or  by  its  voltage  transfer  function* 


^Integrals  and  sums  without  limits  are  over  the  entire  range  (-oo,  -hoo)  of 
nonzero  integrands  and  summands,  respectively. 
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I 


H(f)  =  cfT  exp(-i2irfr)  h(r)  . 


(5) 


Particular  examples  are  the  nonimpulsive  filters, 


h(r)  = 


^  for  0  <  T<  T 


0  otherwise 


,  box  car  filter. 


(6) 


and 


RC 


h(T)  = 


RC  filter 


(7) 


Another  class  of  great  interest  is  the  impulsive  filters  with  a  response 
composed  of  a  number  of  equispaced  impulses: 


h(r) 


w(n)  6(r-nA)  . 


(8) 


Sample  increment  A  is  arbitrary;  the  sum  is  over  all  nonzero  weights  {w(n)]  . 
Since  the  filter  output  in  figure  1  in  steady  state  is  given,  in  this  latter 
case,  by 


z(t)  = 


dr  h(r)  y(t--0  = 


dr  ^  w(n)  6Cr-nA)  g{x(t-T)}  = 


=  ^  w(n)  g{x(t-nA)) 


(9) 


an  equivalent  alternative  form  is  that  given  in  figure  2.  Namely,  the  input 
x(t)  is  sampled  at  increments  A  apart  in  time  and  subjected  to  nonlinearity 
g[x) .  These  quantities  are  then  weighted  and  accumulated,  to  give  output  z(t) 
at  a  time  t  equal  to  a  multiple  of  A.  As  a  special  case,  if 
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w(n)  = 


1  for  1  <  n  <  M 
0  otherwise. 


’  ,  box  car  filter, 


then  all  M  samples  of  the  input  are  equiweighted . 


(10) 


The  input  spectrum  that  will  be  most  closely  considered  is  taken  to  be 
flat  in  a  band  W  about  center  frequency  +f^;  see  figure  3.  The  Q  of  this 
spectrum  is 

0  =  fj-ZW  >  ^  .  (11) 

and  the  highest  frequency  contained  is 

fh  ■  fc  "  2  ■  "  l)  ■ 

The  constraint  on  Q  in  (11)  guarantees  that  W  is  always  the  bandwidth  of  the 
positive  frequency  components  of  the  input  spectrum. 


The  normalized  correlation  corresponding  to  figure  3  is 


p{T)  =  cos(2irf^T)  sinc(Wt)  , 


(13) 


where  we  define 


sinc(x)  = 


sin(-trx) 

irX 


(14) 


For  the  special  case  of  a  low-pass  spectrum,  f^  =  W/2,  (13)  reduces  to 


p{T)  =  sinc(2Wi:)  for  Q  =  -  ; 


(15) 


the  input  spectrum  is  flat  over  (-W,W)  in  this  particular  case, 
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2 

Let  and  denote  the  mean  and  variance,  respectively,  of  the 


output  of  the  detection  systems  in  figures  1  and  2.  Our  main  interest  here  is 
in  the  evaluation  of  the  output  signal-to-noise  ratio  y  defined  in  accordance 
with  the  power  deflection  criterion; 


(16) 


The  subscripts  1  and  0  denote  the  corresponding  hypotheses,  as  already 
introduced  in  (1).  We  will  determine  the  dependence  of  y  on  all  the 
parameters  encountered  above,  such  as  ,  Oq,  PqCH.  v,  T,  a,  M,  W,  0, 
and  compare  the  performance  of  different  nonlinear  systems  for  various 
sampling  rates  and  spectra.  These  results  will  greatly  extend  those  given  in 
[1],  for  example,  and  will  be  much  more  accurate,  since  we  will  use  1000  terms 
in  our  series  expansions,  instead  of  the  3  or  4  terms  used  there.  Observe 
that  the  absolute  levels  (gains)  of  the  nonlinearity  and  filter  cancel  out  in 
quantity  y;  thus,  we  can  assume  any  convenient  level  for  them,  as  done  in  (3), 
(4),  (6),  (10),  for  example.  Limitations  of  the  output  deflection  criterion 
(16)  will  be  discussed  later. 
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DERIVATION  OF  OUTPUT  MOMENTS 


Output  Mean 


In  order  to  determine  the  detection  system  output  signai-to-noise  ratio 
defined  in  (16),  we  need  m^  and  a^.  For  any  filter  h(r),  whether 
impulsive  or  not,  since  the  output  in  steady  state  is 


z(t)  =  dr  h(r)  y(t-T)  , 


then  the  mean  of  the  output  is 


=  z(t)  =  J  dr  h(T)  y(t-T)  =  niy 


r 


dT  h(r)  , 


(17) 


(18) 


in  terms  of  the  mean  of  nonlinearity  output  y(t).  (An  overbar  denotes  an 
ensemble  average.)  However,  this  latter  quantity  is  given  by 


rriy  =  y(t)  =  g[x(t)] 


(19) 


in  terms  of  the  first-order  probability  density  function  p^^^  of  input  x(t). 

2 

Since  the  input  is  zero  mean  Gaussian,  with  variance  a  ,  then 


P 


(1) 


(X) 


where  we  define 


(20) 


0(W) 


(2ir) 


-1/2 


(21) 


exp(-w^/2)  . 
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Substitution  of  (20)  in  (19)  yields  for  the  mean  of  y(t), 


dx  g[x]  -  =  I  dw  g{ow]  0(w) 


(22) 


Coupled  with  (18),  system  output  mean  m^  is  now  available  as 


m. 


=  y  dw  gfow|  ia(w) 


r 


dr  h(r)  . 


Output  Variance 

The  ac  (zero-mean)  output  of  the  detection  system  is 
z(t)  =  z(t)  -  Kt)  =  fdT  h(T)  y(t-T)  , 

V 

where  the  ac  filter  input  is 


y(t)  =  y(t)  -  y(t)  =  y(t)  -  m^  . 

Then  the  variance  of  the  detection  system  output  is 


(23) 


(24) 


=  z^(t)  =  JJ  du  dv  h(u)  h(v)  y(t-u)  y(t-v)  = 

=  JJdu  dv  h(u)  h(v)  R-(u-v)  =  JdTa(r)  R~(r)  . 


(25) 


where 


RyCT)  =  y(t)  y(t-l) 

is  the  (auto)  correlation  of  random  process  7(t),  and 

a(r)  =  J  du  h(u)  h(u-T)  (26) 

is  the  autocorrelation  of  the  deterministic  filter  impulse  response. 
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An  alternative  form  to  (25)  is  available  via  the  frequency  domain 
expression 

=  ^df  |H(f)|  ^  P^(f)  . 


(27) 


where  [H(f)[  ^  is  the  filter  power  transfer  function,  and  Py(f)  "is  the 
power  density  spectrum  of  y(t). 


Now  from  (24),  we  see  that  correlation 

R~(T)  =  R  (T)  -  m^  . 
y  y  y 

Also,  from  (2),  we  have  correlation 


(28) 


Ry(r)  =  y(t)  y(t-r)  =  glx(t)}  glx(t-x)]  = 

=  Jf^^a  9^1  S^^bl  *  (29) 

in  terms  of  the  second-order  probability  density  function  of  input  process 
x(t). 


At  this  point,  we  use  Mehler's  expansion  for  a  Gaussian  process  [2,  (67)]: 


oO 


p(^^(Xa.  Xjj,  D  =  p(^^x^)  p(^^x^) 


2  It 


He 


k\  a , 


(30) 


k=0 


where  p^^^  has  already  been  encountered  in  (20)-(21),  and  He|^(x)  is  the 
Hermite  polynomial  [3,  (22.2.15)].  Substitution  of  (30)  in  (29)  yields 
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Ry(r)  =  ^  p^r)  G^(k)  . 


(31) 


k=0 


where  nonlinearity  coefficients  are  defined  as 


G(k)  = 


dx  g^x]  p^^^(x)  He|^(x/<j)  = 


J 


dw  g[owJ  0(w)  He|^(w)  for  k  >  0  . 


(32) 


Here  we  used  (20)  and  (21). 


In  particular,  we  have  from  (32)  and  (19),  the  zero-th  order  coefficient 


(33) 


2 

Thus  the  k  =  0  term  in  sum  (31)  is  simply  m^.  Combining  this  information 
with  (28)  and  (31),  we  obtain  correlation 


co 


(34) 


k=l 

It  is  worthwhile  observing  that  coefficients  [G(k)]  in  (32)  depend  solely 
on  the  input  standard  deviation  a  and  the  nonlinearity  g  [x]  ;  they  are 
independent  of  the  input  spectral  shape  or  the  filter  characteristics.  Of 
course,  from  (1),  since 


a  under  hypothesis  H 


Oq  under  hypothesis  Hq 


(35) 


J 
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then 


G(l<) 


G^(k)  underH 


'GQ(k)  under  H 


(36) 


However,  we  will  keep  a  general  for  now,  at  least  until  we  have  to  specialize 
to  H-|  versus  Hq,  or  to  evaluate  output  signal-to-noise  ratio  y  in  (16).  A 
similar  procedure  has  been  adopted  with  respect  to  general  p{T)  in  (30),  (31), 
and  (34),  above. 


We  now  utilize  (34)  in  (25)  to  obtain  the  detection  system  output  variance 


=  ^  ^  G^(k)  A(k)  , 

k=l 


(37) 


where  we  define  coefficient 


A(k) 


dr  a(T) 


p\t) 


for  k  >  1  . 


(38) 


This  sequence  [A(k)}  depends  solely  on  filter  h(r)  and  input  normalized 
correlation  p(r) ;  it  is  independent  of  a  and  g^x^.  Thus  sequences  [G(k)]  and 
[A(k)j  in  (32)  and  (38),  respectively,  completely  separate  the  dependence  of 
the  system  output  variance  on  the  relevant  parameters  of  the  problem. 

However,  just  as  in  (1)  and  (36),  since 
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then 


rp^(r)  under 

p(r)  ={  (  . 

Ipq{t)  under' HqJ 


A(l<) 


(k)  under  H 


AQ(k)  under  H 


(39) 


(40) 


The  results  above  for  output  mean  in  (18)  and  output  variance 

2 

in  (37)  are  exact.  There  are  no  assumptions  regarding  small  input 
signal-to-noi se  ratio,  large  averaging  time,  or  large  time-bandwidth  product. 
They  hold  for  arbitrary  input  strength  a,  input  normalized  correlation  p(r) , 
filter  impulse  response  h(T),  and  nonlinearity  glx].  Also  a  and  p(r)  can  vary 
with  the  hypothesis,  as  in  (1).  The  input  spectrum  can  be  low-pass, 
broadband,  or  narrowband.  The  filter  can  be  impulsive  as  in  (8),  or  otherwise 
as  in  (6)  and  (7);  the  sampling  interval  A  and  weights  [w(nj3  in  (8),  and  the 
duration  T  or  time  constant  RC  in  (6)  or  (7),  are  arbitrary.  The  nonlinearity 
g[x]  can  be  a  v-th  law  power  device  as  described  in  (3)  or  (4),  but  need  not 
be;  also,  w  is  not  limited  to  being  an  integer.  The  effects  of  deliberately 
undersampling  the  input  process  in  figure  2  can  be  investigated  by  choosing 
sampling  increment  A  larger  than  the  inverse  of  twice  the  highest  frequency, 
(12),  in  figure  3;  conversely,  the  effects  of  a  continuous  filter  impulse 
response  can  be  deduced  by  choosing  A  very  small.  All  of  these  effects  will 
be  investigated  here. 


The  major  assumption  utilized  is  that  the  input  signal  and  noise  must  be 
Gaussian;  this  precludes  the  presence  of  pure  tones  in  the  input  signal. 


12 


TR  7633 


Output  Siqnal-to-Noise  Ratio 

The  system  output  signal-to-noi se  ratio  was  defined  in  (16).  We  now 
employ  (18),  (33),  and  (37)  to  obtain  it  in  the  form 


(41) 


The  utility  of  this  result  depends  on  the  ability  to  accurately  and 
efficiently  evaluate  the  single  integrals  for  G(k)  and  A(k)  in  (32)  and  (38), 
respectively,  for  high-order  k. 

A  note  of  caution  is  worthwhile  here:  since  output  signal-to-noi se 
ratio,  Y,  in  (16)  and  (41)  only  uses  second-order  moment  information,  it  will 
have  limited  capability  insofar  as  determining  the  system  operating 
characteristics,  that  is,  detection  probability  versus  false  alarm 
probability,  unless  output  z(t)  is  fairly  well  approximated  by  a  Gaussian 
random  variable.  This  latter  situation  will  obtain  when  the  product  of 
averaging  time  and  input  bandwidth  is  large  relative  to  1;  furthermore,  that 
means  that  low  input  signal-to-noise  ratios  can  be  tolerated  and  yet  decent 
performance  predictions  can  be  realized.  Hence,  although  we  concentrate  here 
on  the  statistic  y.  we  are  aware  of  its  limitations  as  a  performance  measure. 
To  accurately  determine  the  exact  operating  characteristics,  the  techniques  of 
[4]  could  be  advantageously  employed. 
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BASIS  OF  COMPARISON 


In  order  to  compare  the  various  nonlinear  systems  subject  to  different 
input  spectra  and  averaging  filters,  a  standard  performance  level  for  output 
signal-to-noise  ratio  y  will  be  adopted.  All  comparisons  will  then  be  made 
with  this  standard,  which  will  now  be  derived. 

The  standard  system  of  interest  is  depicted  in  figure  4;  it  is  a  special 

2 

case  of  figure  1  with  g(x]  =  x  for  all  x,  that  is,  the  full-wave  square-law 
rectifier.  Since  the  output  for  figure  4  is 


(42) 


then  mean 


m^  =  z(t)  =  H(0)  , 


(43) 


where  we  utilized  (5).  The  total  output  power  (ac  and  dc)  is 


X 


(44) 


2 

in  terms  of  the  spectrum  of  process  x  (t).  But  since  the  correlation  of 
2 

X  (t)  IS,  for  Gaussian  process  x(t). 


R  ^(T)  =  x^(t)  x^t-T)  =  +  2o^  p^{T)  , 


(45) 


X 


2 

then  the  spectrum  of  x  (t)  is 


P  ^(f)  =  *(f)  +  P^^^f)  . 


(46) 
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Figure  4.  Full-Wave  Square-Lauj  Detection  System 

f 
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Here,  and  in  the  following,  we  will  use  the  notation 
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=  /  dT:  exp(-i2irfT)  p'^(r) 


(47) 


Thus,  is  the  k-fold  convolution  of 


P^^^f) 


=  J 


dr  exp(-i2irfr)  pix)  • 


(48) 


For  the  flat  bandpass  spectrum  of  major  interest,  P^^^(f)  is  just  the 

2  n  1 

spectrum  of  figure  3,  without  the  factor  a  .  More  generally,  P'  '(f)  is 

(2) 

the  spectrum  of  input  x(t),  normalized  to  unit  area.  P  (f)  is  depicted  in 
figure  5.  Combining  (43),  (44),  and  (46),  we  find  output  variance 


2  o 
=  da 


I  df  |H(f)|  ^  P^^\f)  . 


(49) 


The  output  signal-to-noise  ratio  of  our  standard  system  follows  upon 
using  (43)  and  (49)  in  (16): 


H^(0) 


^s  = 


2o; 


j"  df  |H(f)[ 


(50) 


where  we  have  added  subscripts  to  differentiate  hypotheses  H.|  and  If 

we  wanted  to  maximize  by  choice  of  the  filter,  (50)  indicates  that  the 
relative  power  transfer  function  jH(f)/H(0)j  ^  of  the  low-pass  filter  should 
be  made  as  sharp  (narrow)  as  possible  about  f  =  0,  where  the  power  ratio  must 
necessarily  be  1.  But  since,  in  practice,  the  effective  duration  T  of  filter 


impulse  response  h(T)  is  limited,  there  is  an  upper  limit  to 
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Rather  than  attempt  to  derive  this  absolute  optimum  value  of  we 
develop  an  approximation  to  in  the  case  where  the  filter  width,  1/T,  is 
much  narrower  than  the  input  spectral  width,  W;  that  is,  TW  »  1 ,  which 
corresponds  to  a  long  averaging  time  assumption.  In  this  case,  (50)  yields 
the  approximation 


(51) 


where 


(52) 


is  defined  as  the  effective  duration  of  filter  impulse  response  h(T).  The 
last  quantity  in  (51)  is  the  desired  approximation .  y  ,  to  be  utilized  as  a 
basis  of  comparison.  In  deriving  this  result,  we  have  utilized  (47)  in  the 
form 


(53) 


and  Parseval's  theorem  for  the  Fourier  transform  pair  in  (48). 
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For  the  flat  bandpass  spectrum  of  figure  3,  where  W  measures  the  width  of 
the  positive  frequency  components  of  the  input  spectrum,  that  is, 

p<i)(f)=L  (54) 

the  approximation  in  (51)  becomes 


for  flat  bandpass  spectrum  with  Q 


(55) 


If  the  input  spectrum  is  modified  from  figure  3,  this  result  must  be 
re-evaluated  from  (51).  In  fact,  the  interpretation  of  bandwidth  W  must  be 
done  carefully  and  with  precision. 


In  the  further  special  case  where  hypothesis  Hq  corresponds  to  noise 
only,  and  to  signal  plus  noise,  then 


where  N  and  S  are  the  input  noise  and  signal  powers,  respectively,  and  (55) 
reduces  to 

fsV 

Tb  =  for  flat  bandpass  spectrum  .  (57) 

This,  finally,  is  the  basis  of  comparison,  for  the  standard  system 
s i gna 1 -to-noi se  ratio,  to  be  used  for  all  the  quantitative  results  for  the 
flat  bandpass  spectrum,  and  for  the  two  hypotheses  described  in  (56).  There 
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is  a  certain  amount  of  arbitrariness  in  adopting  (57)  as  a  basis;  however,  if 
two  systems  are  both  compared  with  (57),  then  the  difference  of  those  two 
relative  signal-to-noise  ratios  is  exact,  regardless  of  the  basis. 

It  is  very  important  to  observe  that  the  approximation  for  the  standard 
system  output  signal-to-noise  ratio  in  (51)  is  not  the  optimum  or  maximum 

O 

value  of  the  more  general  result  in  (50)  for  the  square-law  system 
depicted  in  figure  4.  Rather,  Yc  will  be  greater  than  y,  in  some  cases. 

To  see  this,  the  terms  in  the  system  output  variance  (denominator  of  (50)), 
are  illustrated  in  figure  5.  It  is  obvious,  since  everything  is  nonnegative, 
that 


(58) 


Since  this  right-hand  side  is  just  the  quantity  in  the  denominator  of  the  top 
line  of  approximation  (51),  this  means  that  y  >  T  • 

S  d 

More  generally,  this  means  that  we  have  to  expect  the  possibility  that 
the  general  system  output  signal-to-noise  ratio  in  (16)  will  have  y  >  in 
some  cases.  Nevertheless,  because  of  its  simpl icity,  Yj^  in  (57)  will  be  kept 
as  our  basis  of  comparison  for  the  various  systems.  However,  the  interpreta¬ 
tion  of  T  and  W  must  be  carefully  noted  for  each  case. 

The  exact  way  in  which  we  use  the  basis  y^  in  (57)  is  as  follows:  for 
the  general  nonlinear  system  in  figure  1,  we  set  its  output  signal-to-noise 
ratio  Y  equal  to  the  basis,  that  is,  set 
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Y  = 


(59) 


where  S/N  is  the  actual  input  signal-to-noise  ratio  to  the  general  nonlinear 
system,  and  (S/N)^  is  the  input  signal-to-noise  ratio  to  the  standard 
system.  Then  we  solve  (59)  for  the  required  input  signal-to-noise  ratio  S/N 
to  achieve  this  level  of  performance,  and  compute  the  decibel  difference  at 
the  input. 


dB  =  10  log 


S/N 

(S/N) 


(60) 


relative  to  the  standard  system  input  signal-to-noise  ratio. 


In  general,  this  quantity  will  be  a  function  of  the  standard  system  input 
signal-to-noise  ratio  (S/N)^;  however,  for  low  input  signal-to-noise  ratio, 
it  is  independent  of  (S/N)^.  It  affords  a  measure  of  how  much  more  input 
signal-to-noise  ratio  is  required  for  the  general  nonlinear  system  of 
interest,  relative  to  the  square-law  standard  of  figure  4.  In  keeping  with 
the  discussion  above,  we  can  expect  that  the  quantity,  dB,  in  (60)  will  become 
negative  in  some  cases;  this  simply  means  that  the  performance  of  that 
particular  nonlinear  system  is  somewhat  better  than  the  arbitrary  basis 
adopted  in  (57). 
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OUTPUT  SIGNAL-TO-NOISE  RATIO  FOR  TWO  FILTER  CLASSES 

The  system  output  signal-to-noise  ratio  y  was  derived  in  (41)  in  terms  of 
the  general  filter  impulse  response  h(X).  In  this  section,  we  shall  develop  y 
in  more  detail  for  impulsive  filters,  as  in  (8),  and  nonimpulsive  filters,  as 
in  (6)  and  (7).  The  latter  case  correponds  to  an  analog  filtering  procedure. 

Impulsive  Filters 


The  impulse  response  h(T)  takes  the  form  (8)  in  this  case,  and  the 
corresponding  system  block  diagram  is  given  in  figure  2.  We  substitute 
(8)  in  (26)  to  obtain  the  autocorrelation  of  h(T)  as 

a(T)  =  ^  b(n)  6(T-nA)  ,  (61) 

n 

where  sequence  ^b(nj}  is  the  autocorrelation  of  the  filter  weights; 

b(n)  =  ^  w(m)  w(m-n)  for  all  n  .  (62) 

m 

The  evaluation  of  A(k)  in  (38)  is  immediate  in  this  case,  yielding 

A(k)  =  I  n)  p*^(nA)  for  k  >  1  .  (63) 

n 

There  follows  from  (18),  (33),  and  (8),  the  system  output  mean  in  the  form 

m^  =  G(0)  ^w(n)  ,  (64) 

n 
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while  from  (37)  and  (63),  the  variance  is  given  by 


Od 


It  G^(k)  2  b(n)  p^nA)  . 


(65) 


k=l 


Combining  these  last  two  expressions  in  (16),  the  output  signal-to-noise 
ratio  for  the  impulsive  class  of  filters  can  be  expressed  as 


Yj  = 


M  [g^(0)  -  Gq(0)]' 


^  G^(k)  ^  B(n)  p^(nA) 


(66) 


k=l 


where 


M  = 


^  w(n)'] 
.  n  J 


^w^  ( n ) 


(67) 


is  the  effective  number  of  samples  in  impulsive  filter  response  h(T)  in  (8), 
and 


w(m)  w(m-n) 


^  kinl  ^  _iL 

^  ^  ^  w2(m) 


(68) 


m 


is  the  normalized  autocorrelation  of  the  filter  sampling  weights  [w(n)]  .  The 
result  (66)  holds  for  any  impulsive  filter  with  any  weight  sequence  [w(n))  , 
sampling  increment  A,  input  levels  ,  Oq,  input  correlation  PQ(Tr).  and 
nonlinearity  g {x]  . 
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Nonimpulsi ve  Filters 


Examples  of  filter  impulse  responses  in  this  class  were  presented  in 
(6)  and  (7).  In  the  former  case,  the  autocorrelation  (26)  is 


a(T)  = 


}  (l  -  for  |t1<  T 


0  otherwise 
while  in  the  latter  case, 

1  f  111 


box  car  filter. 


(69) 


a(r)  =  j 


RC 


exp 


(-1) 


for  al 1  T  ,  RC  filter 


(70) 


If  we  substitute  (38)  in  (41),  the  output  signal-to-noise  ratio  for  the 
nonimpulsive  class  of  filters  can  be  expressed  as 


T 


[g^(0)  -  Gq(0)] 


7'w  TO  •)  2  r  k 

GQ(k)  dTa(T)  Pq(T) 

k=l 


(71) 


where 


T  = 


f' 

fdr  h(r)] 

J 

[dr  h^(T) 

(72) 


is  the  effective  duration  of  filter  impulse  response  h(r),  and 

J  du  h(u)  h(u-T) 

^  '  f,  ,2,  ' - 

'  J  du  h  ( u) 

is  the  normalized  autocorrelation  of  the  filter  response.  The  result 
(71)  holds  for  any  input  levels  a^,  Oq,  input  correlation  pQ(r) ,  and 
nonlinearity  g(x].  However,  it  does  not  cover  a  filter  containing  any 


(73) 
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impulses,  since  the  denominator  integral  in  (72)  and  (73)  is  infinite  then. 
Nevertheless,  the  corresponding  result  in  (66)  for  the  impulsive  filters  can 
be  derived  as  a  limit  of  a  set  of  progressively  narrower  pulses  at  multiples 
of  increment  a.  Thus,  (71)  does  have  the  capability  of  covering  the  most 
general  filter  structure,  if  manipulated  properly. 

Long  Averaging  Time  Assumption 


The  results  in  (66)  and  (71)  for  the  output  signal-to-noi se  ratio 
simplify  somewhat  when  the  averaging  time  of  the  filter  (Ma  or  T)  is  much 
larger  than  the  correlation  time  (1/W)  of  the  input  process  x(t)  under 
hypothesis  Hq.  In  the  impulsive  filter  case,  this  means  that  the  sum  on  n 
in  the  denominator  of  (66)  can  be  approximated  according  to 


^  B(n)  pQ(na)  s  B(0)  ^  p^(na)=^  PQ(na)  for  k  >  1  .  (74) 

n  n  n 


leading  to  approximation 


M  [g^(0)  -  Gq(0)]^ 


(75) 


It  is  interesting  to  observe  that  the  exact  detailed  values  of  the  weights 
[w(n)]  in  impulsive  filter  response  (8)  are  immaterial  to  the  value  of  (75), 
except  insofar  as  they  affect  effective  number  M  via  (67).  This  simplifi¬ 
cation  is  not  possible  for  the  general  averaging  time  result  in  (66),  which 
depends  on  the  weights  through  their  normalized  autocorrelation  sequence 
[B(nj^  in  (68). 
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On  the  other  hand,  for  nonimpulsive  filters,  the  long  averaging  time 
assumption  means  that  the  integral  on  T  in  the  denominator  of  (71)  can  be 
approximated  according  to 


J dr  a(r)  -  “(0)  J  dr  pqCD  =  J  cit  p^(r) 


for  k  >  1 


(76) 


leading  to  approximation 


T  [g,(0)  -  Gg(0)72 

.sa.  1  2  r  k  ■ 

k=l 


(77) 


Analogous  to  the  observation  above,  the  shape  of  the  detailed  impulse  response 
h(r)  is  immaterial  to  the  value  of  (77),  except  as  it  affects  the  effective 
duration  T  via  (72) . 


The  approximations  in  (75)  and  (77)  for  the  long  averaging  time 
assumption  are  not  used  in  the  numerical  results  that  follow  later.  Rather, 
the  exact  result  (66)  for  impulsive  filters, and  (71)  for  nonimpulsive  filters, 
are  extensively  utilized.  Also  the  danger  of  using  the  long  averaging  time 
assumption  when  inappropriate  is  illustrated  by  a  numerical  example  in  a  later 
section. 


Narrowband  Input  Spectrum 

We  now  return  to  general  output  signal-to-noise  ratio  y,,  for  arbitrary 

N 

averaging  time  in  (71)  and  consider,  for  the  moment,  the  case  of  a  very 
narrowband  input  spectrum  under  Hq.  That  is,  let  normalized  correlation 
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Po(T)  =  r(T)  cos(2^f^T+  (J(T))  ,  (78) 

where  center  frequency  is  much  larger  than  the  highest  frequency  contents 
of  envelope  r(T)  and  phase  0(t);  this  is  called  a  high-Q  input  spectrum.  Then 
the  integral  in  the  denominator  of  (71)  becomes 


dr  a(r)  pq(t)  =  Jdra(T)  r^r)  cos^2irf^r+  0(t))  = 


k 

cos  © 


0  for  k  odd 


dr  a(r)  r'"(T)  = 


2^  Vk/2/  J 


dT  a(r)  r^(r)  for  k  even 


(79) 


j 


Here,  the  k-th  power  of  the  cosine  varies  so  quickly  with  T  that  we  replaced 
it  by  its  average  value  and  removed  it  from  under  the  integral  on  T. 
Substitution  of  (79)  in  (71)  yields  the  system  output  signal-to-noise  ratio 


T  [g,(0)  -  Gp(0)]' 


'NB 


ooZ 


FT  ik 


k=2 

even 


drcx(T)  r^T) 


for  a  narrowband  input  spectrum. 


(80) 


This  form  for  y,,_  leads  to  an  interesting  conclusion  regarding 
Nd 

symmetric  full-wave  rectifiers  versus  half-wave  rectifiers.  Namely,  reference 
to  the  defining  relation  (32)  for  nonlinearity  coefficient  G(k)  reveals  that, 
for  k  even,  a  symmetric  full-wave  rectifier  has  a  value  of  G(k)  exactly  double 
that  for  the  corresponding  half-wave  rectifier.  However,  since  (80)  only 
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involves  the  even  k  values,  whether  a  full-wave  rectifier  or  half-wave 
rectifier,  this  factor  of  2  cancels,  and  output  signal-to-noise  ratio  in 

(80)  is  exactly  the  same  for  a  symmetric  full-wave  rectifier  as  for  a 
half-wave  rectifier.  This  conclusion  holds  for  any  rectifier  g{x],  not  just 
the  u-th  law  rectifiers  in  (3)  and  (4).  It  also  holds  for  any  filter  h(r)  and 
normalized  correlation  envelope  and  phase  r(r)  and  0(T)  in  (78),  and  is  not 
limited  to  large  averaging  times  or  small  input  signal-to-noise  ratio.  The 
only  restriction  is  the  required  high  Q  of  the  input  spectrum. 

This  conclusion  regarding  identical  values  for  symmetric  full-wave 
rectifiers  and  half-wave  rectifiers  is  also  physically  reasonable  for  a  high-Q 
input,  in  that  no  significantly  different  information  is  contained  by  the 
negative  lobes  of  the  waveform  x(t)  when  its  envelope  and  phase  are  slowly 
varying  relative  to  the  center  frequency. 
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w-TH  LAW  RECTIFIERS 

Up  to  this  point,  all  the  relations  involving  nonlinearity  coefficient 
G(k),  defined  in  (32),  have  been  general.  We  now  specialize  to  the  symmetric 
full-wave  and  half-wave  w-th  law  rectifiers  mentioned  in  (3)  and  (4), 

Full-Wave  Rectifiers 

The  nonlinearity  of  interest  here  is 

g{x}  =  )xl "  for  all  x  (v  >  0)  ,  (81) 

which  is  a  symmetric  full-wave  rectifier.  Substitution  of  (81)  in  (32)  yields 

G(k)  =  ^  Iwl t5(w)  He|^(w)  for  k  >  0  ,  (82) 

where  0(w)  is  defined  in  (21),  and  He|^(w)  is  a  Hermite  polynomial 

[3,  (22.2.15)].  Since  He|^(w)  is  odd  in  w  for  k  odd,  it  immediately  follows 

that 

G(k)  =  0  for  k  odd  .  (83) 

On  the  other  hand, 

OO 

G(k)  =  2a^  J  dw  w’^  0(w)  He|^(w)  for  k  even  (0,  2,  4,  ...)  .  (84) 

0 
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In  order  to  evaluate  this  integral,  we  first  define  the  quantity 

oo 

L(k)  =  j'  dw  w"'  0(w)  He|^(w)  for  all  k  >  0  .  (85) 

0 

It  should  be  noted  that  L(k)  is  a  function  of  v,  in  addition  to  the  explicitly 
indicated  dependence  on  k;  however,  it  is  independent  of  a.  A  recurrence  for 
L(k)  is  derived  in  appendix  A: 


L(k)  =  (v+2-k)  L(k-2)  for  all  k  >  2  , 


(86) 


with  starting  values 


^-1 


L(0)  =2^  ir  ^  ,  L(l)  =  2  ^  ir 


(87) 


Thus  very  high  order  values  of  L(k)  can  be  quickly  evaluated  with  the  aid  of 
just  two  gamma  function  computations. 

The  nonlinearity  coefficient  in  (84)  is  therefore  expressible  as 


or 


G(k)  =  2a  L(k)  for  k  =  0,  2,  4 . 


G(k)  =  (v+2-k)  G(k-2)  for  k  =  2,  4,  6, 


(88) 


(89) 


with  starting  value 
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We  now  utilize  the  results  in  (88)  and  (83)  to  express  the  impulsive fi Iter 
system  output  signal-to-noise  ratio  in  (66)  as 


^  U(k)  ^  B(n)  p^(na) 

fe  n 

even 


(91) 


where  subscripts  IF  denote  the  impulsive  filter  full-wave  rectifier  case  and 
where  we  have  also  defined  sequence 


U(k)  =  Jt  fo*"  an  k  >  0  .  (92) 

''•  L'^(O) 

This  latter  sequence  has  a  simple  recurrence,  as  seen  by  reference  to  (86), 
namely 


U(k)  =  U(k-2)  for  all  k  >  2  , 


(93) 


with  starting  values,  from  (87), 


U(0)  =  1  ,  U(l)  =  2 


(94) 


In  fact,  since  (91)  only  involves  even  k,  can  actually  be  evaluated 
without  the  aid  of  any  gamma  functions. 


An  exactly  analogous  procedure  applied  to  nonimpulsive  result  (71)  leads 
to  output  signal-to-noise  ratio 
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'NF 


T 

(.<'<  - 

oc 

2 

m)  ]■ 

dra(T)  Po(T) 

k=2 

even 


(95) 


where  subscripts  NF  denote  a  nonimpulsive  filter  and  a  symmetric  full-wave 
rectifier. 


Full-Wave  Square-Law  Rectifier 


Here  we  specialize  the  above  results  to  the  case  of  u  =  2,  that  is, 
2 

=  X  for  all  x.  There  follows,  upon  use  of  (93)  and  (94),  the 
significantly  simpler  results 


T,p(»  -  2) 


.  "  (A'%  - ^ 

2  ^  B(n)  p2(nA) 


from  (91 ) ,  and 


^NF^"  -  2) 


[dr. 


(1^)  Po(r) 


(96) 


(97) 


from  (95).  This  fortuitous  situation  occurs  because  the  recurrence  (93) 
generates  zero  coefficients  for  k  >  4,  when  «  =  2.  More  generally,  v  even 
would  also  terminate  the  recurrence  at  k  =  v  +  2. 
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Full-Wave  Linear  Rectifier 


This  special  case  corresponds  to  v  =  1 ,  that  is,  gfx]  =  1x1  for  all  x. 


Then 


T  j  p  ( 1  )  oo. 


(98) 


U(k)  >  B(n)  pJCna) 

k=2  n 

even 


from  (91),  where  fu(k)}  is  given  by  (93)  with  u  =  1.  Also,  from  (95), 


Ynf(''  =  1)  = 


OQ 

k=2 

even 


h 


u(k)  J  dr  a(r)  pq(t:) 


Half-Wave  Rectifiers 

We  now  return  to  general  values  of  power  law  w,  but  to  half-wave 
rectifiers  characterized  by 


g[x}  = 


X  for  X  >  0 


0  for  X  <  0 


Substitution  of  (100)  in  (32)  yields  nonlinearity  coefficient 


(99) 


(100) 


G{k)  =  o  L(k)  for  all  k  >  0  , 


(101) 
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where  L(k)  was  defined  in  (85).  When  (101)  is  utilized  in  (66),  the  result 
for  the  output  signal-to-noise  ratio  is 


(102) 


k=l 


n 


where  IH  denotes  impulsive  filters  and  half-wave  rectifiers.  This  result  is 
identical  to  Yjp  in  (91),  except  for  the  inclusion  of  all  the  odd  k  values 
here.  Similarly,  (71)  leads  to 


(103) 


for  nonimpulsive  filters  and  half-wave  rectifiers,  NH.  This  is  identical  to 

in  (95)  except  for  the  inclusion  of  all  the  odd  k  values  here,  in  the 
Nr 

case  of  half-wave  rectification. 

Half-Wave  Square-Law  Rectifier 

When  we  set  u  =  2  in  the  results  above,  and  refer  to  (93),  we  find  that 


Yih(«  =  2)  = 


"  (°^°o  - 


2 


(104) 


k=l 

odd 


n 


n 
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while 


=  2)  = 


O  oO 


2  J  dr  a(-r)  PQ(r)  +  ^  U(k)  J  dr  a(r)  pQ(r) 

k=l 

odd 


(105) 


Neither  of  the  series  in  odd  k  terminate,  since  the  numerator  in  recurrence 
2 

(93)  is  (4-k)  ,  which  is  always  nonzero  for  k  odd. 


Half-Wave  Linear  Rectifier 


Upon  setting  u  =  1  in  (93)  and  (94),  we  obtain 


U(0)  =  1  ,  U(l)  =  I  ,  U(k)  =  U(k-2) 


for  k  =  2,  4,  6, 


(106) 


with  all  other  U(k)  zero.  Then  (102)  and  (103)  yield 


Tih(«  =  1) 


OO 


I  ^  B(n)  pQ(nA)  +  ^  U(k)  ^  S(n)  pQ(nA) 

n  k=2  n 

even 


(107) 


and 


^  /  .  -  j)' 


iJdra(T)  Pg(r)  +  ^  U(k)  drcc(r)  p^(r) 


k=2 

even 


(108) 


respecti vely . 
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The  programs  to  be  furnished  later  are  not  limited  to  the  special  cases 
in  (96)-(99)  and  (104)-(108),  but,  in  fact,  cover  arbitrary  values  of  v. 


Equality  of  Performance  for  v  =  1 


It  is  interesting  to  observe  that  these  last  two  results  for  the  linear 
half-wave  rectifier  are  identical  to  the  corresponding  earlier  results  for  the 
symmetric  linear  full-wave  rectifier  in  (98)  and  (99),  except  for  the 
additional  term  for  k  =  1  here.  Thus,  for  nonimpulsive  filters, 


if 


I 


dT  a(r)  PQ(r) 


0 


(109) 


in  (108),  then  the  output  signal-to-noise  ratios  of  the  linear  symmetric 
full-wave  rectifier  and  linear  half-wave  rectifier  are  the  same.  But  since, 
from  appendix  B, 


dra(r)  Pn(r-) 


J 

1  df 

|H(f: 

]■ 

df  |t 

i(f) 

(110) 


the  only  way  (109)  can  be  true  is  if  the  filter  power  transfer  function 
jH(f)[‘^  and  the  input  spectrum  under  Hq,  namely  Pq  '(f),  do  not 
overlap;  see  figure  6  for  the  flat  bandpass  input  spectrum  example.  Then 
approximately,  if 


i .e. ,  0 


(ill) 


then  (109)  is  substantially  satisfied  and  the  linear  half-wave  rectifier  and 
the  symmetric  linear  full-wave  rectifier  have  similar  output  signal-to-noise 
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Figure  6.  Spectral  Quantities  -Tor  (110)  and  (112) 


ratios.  This  lower  bound  on  Q  in  (111),  relative  to  the  inverse  time- 
bandwidth  product,  will  be  encountered  again  later,  when  we  investigate 
general  v-th  law  rectifiers  for  general  bandpass  spectra,  as  delineating  a 
distinct  dichotomy  in  performance  for  half-wave  rectifiers. 

It  will  be  recalled,  in  an  earlier  section  dealing  with  narrowband  input 
spectra,  (78)-(80),  that  the  output  signal-to-noise  ratios  of  a  u-th  law 
half-wave  rectifier  and  a  symmetric  full-wave  rectifier  were  identical  for  all 
w,  if  the  0  of  the  input  spectrum  was  very  large.  The  result  here, 
particularly  (111),  is  a  much  milder  requirement  that  achieves  the  same 
result,  but  only  for  the  linear  device,  v  =  1.  The  reason  that  a  much  more 
moderate  requirement  on  Q  will  suffice  for  v  =  1  is  that  all  the  odd  terms  for 
k  >  1  are  already  absent  from  (107)  and  (108),  whereas  they  were  suppressed  in 
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(79)  only  by  the  use  of  the  very  large  Q  assumption.  So,  for  a  linear 
rectifier,  virtual  equality  of  performance  for  a  half-wave  rectifier  and  a 
symmetric  full-wave  rectifier  can  be  expected  for  a  very  moderate  value  of  Q. 
Larger  values  of  v  would  require  larger  Q  values  for  the  same  result  to 
obtain.  These  observations  will  be  borne  out  by  numerical  results  to  follow 
later. 


For  the  impulsive  filter  result  of  (107),  the  analogous  result  to  (110) 
is  also  derived  in  appendix  B,  giving 


A  ^B(n)  PQ(nA) 
n 


jdf 

Jdf 


where  the  filter  voltage  transfer  function  is  now 


H(f) 


=  2 


w(n)  exp(-i2irfAn) 


(112) 


(113) 


from  (5)  and  (8).  The  integral  in  the  denominator  of  (112)  is  over  any 
frequency  interval  of  length  1/a,  which  is  the  period  in  f  of  the  filter 
voltage  transfer  function  H(f).  The  numerator  integral  in  (112)  is  over  all 
f,  but  could  be  limited,  if  desired,  to  the  fundamental  frequency  interval 
(-1/2A,  1/2A)  if  normalized  input  spectrum  PQ^^(f)  is  replaced  by  its 
aliased  version 


n 


where  the  sum  is  over  all  n  from  -  “o  to  +ao  . 


(114) 
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In  any  event,  (112)  will  be  zero  if  the  filter  power  transfer  function 
2  Cl) 

|H(f))  and  input  spectrum  Pg  ^(f)  do  not  overlap.  Figure  6  is  again 

I  2 

applicable,  but  now  with  the  addition  of  major  lobes  of  |H(f)|  at  f  =  n/a 
for  all  n.  Not  only  must  (111)  be  met  in  order  to  avoid  overlap  near  f  =  0, 
but  the  other  major  lobes  of  jH(f)|  ^  must  not  overlap  PQ^^(f).  This 
requirement  can  be  met,  for  example,  by  keeping  the  lower  skirt,  of  the  lobe 
of  jH(f)|  ^  centered  at  f  =  1/a,  above  the  highest  frequency  of  the  input 
spectrum: 

that  is,  sampling  frequency  fj  =  1/a  must  satisfy 


W 


> 


^  +  0  +  1 

TW  2 


(116) 


This  requirement  is  in  addition  to  that  of  (111). 


Actually,  (116)  is  a  sufficient  condition  for  non-overlap  in  figure  6, 

but  is  not  always  necessary.  In  particular,  for  larger  Q,  the  possibility 

exists  of  deliberately  undersampling  (violating  (115))  and  yet  achieving 

nonoverlap  of  the  aliased  spectral  components  of  (114)  with  the  filter  lobe  at 

f  =  0.  This  will  be  demonstrated  by  example  in  the  results  section,  not  only 

1 

for  V  =  1,  but  for  other  «  values  as  well.  Here  we  have  resorted  to  the 
equality 


Jdf  |H(f)|^  pjj^^f)  =  J^df  |H(f)|  ^  P^^^f)  , 


(117) 


'/4 

based  upon  the  use  of  aliased  spectrum  (114)  and  the  periodic  character  of 
filter  H(f). 
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INDEPENDENT  SAMPLES 


This  section  will  address  the  discrete  detection  system  of  figure  2, 
where  sampling  increment  A  in  (8)  is  taken  such  that  the  samples  of  Gaussian 
process  x(t)  are  statistically  independent.  This  can  be  accomplished  by 
taking  a  large,  in  general.  However,  at  least  in  the  special  case  of  a  flat 
low-pass  spectrum  under  Hq,  as  in  figure  3  with 

=  ^  .  Q  =  ^  ^  for  ifl  <  W  ,  (118) 

then  the  choice  A  =  (2W)  ^  in  corresponding  normalized  correlation  (15)  gives 

[l  for  n  =  o) 

PQ(nA)  =  sinc(n)  =  <  |,  (119) 

|0  otherwise) 

which  also  means  independent  samples.  (Another  special  case  is  afforded  in 
the  flat  bandpass  spectrum  case,  by  choosing  A  =  irrespective  of  center 

frequency  f^;  see  (13)  and  refer  to  (119).) 


For  this  section  only,  dealing  exclusively  with  independent  samples,  we 
will  not  yet  specialize  to  the  case  of  v-th  law  rectifiers,  but  temporarily 
allow  general  nonlinearities.  Then  general  result  (66)  for  the  system  output 
signal-to-noise  ratio  reduces  to 


Tt  = 


M  [g^(0)  -  Gq(0)] 


1 

k=l 


1 


kl  «0<k> 


for  independent  samples. 


(120) 
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by  use  of  (119).  We  observe  immediately  that  the  only  way  the  filter  weights 
{w(nj|  enter  this  result  is  via  effective  number  M  defined  in  (67),  whether 
large  or  small;  normalized  autocorrelation  sequence  {B(n)j  in  (68)  is  not 
relevant  in  this  special  case. 


Numerical  evaluation  of  the  particular  infinite  series  in  the  denominator 
of  (120)  can  be  circumvented,  as  follows;  from  (34)  and  (28), 

oo 

R^iO)  =  ^  ^  G^k)  =  cjy  .  (121) 

k=l 

Then,  also  using  (33),  (120)  becomes 


"  ("vl  -  "'vo)^ 
■  2 
‘’yo 


for  independent  samples  . 


(122) 


But  these  latter  quantities  can  be  obtained  directly  from  the  moments  of  the 
nonlinearity  output. 


niy  =  y(t)  ="g{x(t)]  =  J  dx  g{x}  p^^\x)  , 

y^(t)  =  g^[x(t)]  =  J dx  g^{x]  p^^^x)  ,  (123) 

by  numerical  integration  if  necessary,  once  the  nonlinearity  g[x]  and 
first-order  input  probability  density  function  p^^^(x)  are  specified. 
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Symmetric  Full-Wave  Rectifier 

We  now  specialize  to  symmetric  w-th  law  full-wave  rectifiers  as  in  (4) 
and  a  Gaussian  input  as  in  (20).  (123)  immediately  yields  moments 


„  V  _v/2  -1/2 

m^  =  a  2  ir 


2...  2v  _u  -1/2 
y  (t)  =  o  2  ir 


r(f-) , 


(124) 


Appropriate  substitution  into  (122)  gives  output  signal-to-noise  ratio 


" 

'if  D(w)  -  1 


for  independent  samples  , 


(125) 


where  we  define  quantity 


D(v)  = 


1/2 


(126) 


2  J 


Equation  (125)  is  a  compact  expression  for  the  system  output  signal-to-noise 
ratio  for  any  u  >  0,  requiring  no  summations. 

Special  cases  of  (126)  are 


0(1)  =  I  .  0(2)  =  3  .  0(3)  =  ^  .  0(4)  =  ^ 


(127) 


Thus,  for  example,  (125)  yields 


Tip(v  2)  2 


2 


(128) 
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in  agreement  with  (96)  in  this  case  of  independent  samples. 

More  generally,  since  from  (91)  and  (119),  we  have  output  signal-to-noise 

ratio 


« 

Ytp  =  — ^ —  for  independent  samples  ,  (129) 

2  U(k) 
k=2 
even 

then  comparison  with  (125)  yields  the  following  identity  on  the  [u(k)^ 
sequence  in  (92)-(94): 


<30 

2  udc) 

k=0 

even 


(130) 


This  result  has  been  confirmed  directly,  by  letting  m  =  k/2,  converting  the 
recursion  on  {u(k)}  to  one  on  {U(2m)J  ,  recognizing  it  as  a  hypergeometric 
function,  and  using  [3,  (15.1.20)].  This  identity  will  prove  to  be  very 
useful  in  the  accurate  numerical  evaluation  of  general  result  (91)  for 
symmetric  full-wave  rectifiers  with  statistically  dependent  samples. 


We  now  consider  the  case  where  hypothesis  Hq  corresponds  to  noise-only 
and  to  signal  plus  noise,  that  is. 


2 

a 


=  S  +  N  under  H, 


N  under  H, 


(131) 
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where  N  and  S  are  the  independent  input  noise  and  signal  powers, 
respectively.  When  utilized  in  (125),  we  have  output  signal-to-noise  ratio 


which  is  valid  for  independent  samples  and  all  values  of  M,  v,  and  S/N. 

The  basis  of  comparison  (in  this  section  alone,  for  independent  samples) 
will  be  the  full-wave  square-law  detector,  which  simply  corresponds  to  taking 
V  =  2  in  (132): 


^b 


(133) 


In  keeping  with  the  discussion  surrounding  (59),  we  now  equate  these  last  two 
expressions  above: 


M 


D(v)  -  1 


1  + 


sy/2 

n; 


12 


-  1 


M  fS 
2  N, 


(134) 


where  (S/N)^  is  the  input  signal-to-noise  ratio  to  the  square-law  system, 
and  S/N  is  that  for  the  general  w-th  law  symmetric  full-wave  rectifier.  The 
solution  of  (134)  is 


(135) 
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This  is 
achieve 
system, 
effecti 


the  input  signal-to-noise  ratio 
the  same  output  signal-to-noise 
for  the  independent  samples  case 
ve  number  of  samples. 


requi red 
ratio  as 
.  It  is 


for  the  v-th  law  system  to 
for  the  full-wave  square-law 
independent  of  M,  the 


Although  generally  a  function  of  the  actual  square-law  input 
signal-to-noise  ratio  (S/N)^,  in  the  case  of  low  input  signal-to-noise 
ratio,  (135)  becomes  approximately 

I  a  1  (20(.)  -  2)'^^  for(|]^«l.  (136) 

This  leads  to  the  factor 

Fp(v)  =  I  (20(v)  -  2)^^^  (137) 

by  which  the  input  signal-to-noise  ratio  must  be  increased  for  the  w-th  law 
device  relative  to  the  full-wave  square-law  system,  and  which  is  independent 
of  the  actual  input  signal-to-noise  ratio  (S/N)^.  Subscript  F  denotes 
symmetric  full-wave  rectifiers. 


u 


The  factor 
in  figure  7. 


Fp(v)  is  tabulated  in 
They  both  reveal  that 


decibels  in  table  1  and  plotted  versus 
for  low  input  signal-to-noise  ratio. 
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Table  1.  Additional  Input  Signal-to-Noise  Ratio  Relative  to  Full-Wave 
Square-Law  System;  Low  Input  Signal -to-Noi se  Ratio, 

Symmetric  Full-Wave  Rectifiers,  Selected  v.  Independent  Samples 


dB,  Relative  SNR  dB,  Relative  SNR 
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Figure  7.  Additional  Input  Si gnal -to-No i se 
Ratio  Relative  to  Square-Law  System;  Symmetric 
FWRs ,  Low  Input  Si gnal -to-No i se  Ratio 


Figure  8.  Additional  Input  S i gn a  1 -to-No i se 
Ratio  Relative  to  Square-Law  System; 
Symmetric  Full-Wave  Rectifiers 
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the  fuH-wave  square-law  detector  is  best.  The  linear  full-wave  rectifier 
(FWR),  for  example,  requires  .29' dB  more  input  signal-to-noise  ratio  than  the 
square-law  full-wave  rectifier,  in  order  to  achieve  the  same  output 
signal-to-noise  ratio  (SNR). 

More  generally,  exact  result  (135)  is  plotted  in  figure  8,  for  various 
values  of  the  full-wave  square-law  input  signal-to-noise  ratio  (S/N)^.  The 
curve  labelled  (S/N)^  0  is  exactly  that  already  plotted  in  figure  7.  The 

other  results  for  larger  (S/N)^  seem  to  indicate  that  better  performance  can 
be  achieved  by  choosing  v  larger  than  2,  for  example,  v  =  3.5  for  (S/N)^  = 

0.5.  However,  this  conclusion  is  completely  spurious,  because  the  square-law 
detector  is,  in  fact,  the  optimum  device  to  use  in  this  particular  case  of 
independent  samples  of  the  input;  furthermore,  the  input  samples  should  all  be 
equi-weighted .  These  conclusions  are  based  on  the  likelihood  ratio  derivation 
in  appendix  C.  This  situation  serves  to  accentuate  the  earlier  caution  that 
the  output  signal-to-noise  ratio  of  a  system  with  nonlinearities  is  not  a 
complete  descriptor  of  performance,  and  care  must  be  exercised  in  the  use  and 
interpretation  of  the  system  output  signal-to-noise  ratio;  see  (41)  et  seq. 

The  results  in  figure  8  for  the  higher  signal-to-noise  ratios  cannot  be  used. 

Half-Wave  Rectifier 


We  now  return  to  general  results  (122)  and  (123),  and  consider  half-wave 
rectifiers  as  described  in  (3)  and  a  Gaussian  input  as  in  (20).  Substitution 
in  (123)  yields  moments  which  are  half  of  the  values  listed  in  (124).  When 
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employed  in  (122),  there  follows  for  the  output  signal-to-noise  ratio, 

2 

'IH  -  independent  samples  ,  (138) 


M 


2D(v)  -  1 


where  D(w)  was  defined  in  (126).  For  example, 

VO 


(139) 


for  the  square-law  half-wave  rectifier;  compare  this  result  with  (128)  for  the 
full-wave  rectifier. 


Generally,  from  (102)  and  (119),  we  have 


1 

k=l 


for  independent  samples 


(140) 


Comparison  with  (138)  yields  the  identity 


oa 


2;  U(k)  =  2D(v)  . 
k=0 

Coupled  with  (130),  this  yields 


(141) 


oO 


U(k)  =  D(v)  = 


k=l 

odd 


Twy 


(142) 


Thus,  both  the  even  and  odd  sums  on  give  exactly  the  same  value,  D(w). 

The  result  in  (142)  has  been  confirmed  directly  by  letting  m  =  (k-l)/2, 
converting  the  recursion  on  [U(k)3  to  one  on  fu(2m+l)},  recognizing  it  as  a 
hypergeometric  function,  and  using  [3,  (15.1.20)].  This  identity  will  prove 
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to  be  very  useful  in  the  accurate  numerical  evaluation  of  general  result  (102) 
for  half-wave  rectifiers  with  statistically  dependent  samples. 


We  now  again  consider  the  case  where  hypothesis  corresponds  to 
noise-only  and  to  signal  plus  noise;  see  (131).  Then  (138)  for  the 
output  signal-to-noise  ratio  becomes 


_  M 

^IH  “  2D(w)  -  1 


(143) 


which  is  valid  for  independent  samples  and  all  values  of  M,  v,  and  S/N. 


For  small  input  signal-to-noise  ratio,  this  is  approximately 


=  M 

^IH  ■  2D(v)  -  1 


(144) 


When  this  is  compared  with  the  basis  in  (133)  for  the  ful 1 -wave  square -law 
detector,  we  see  that  the  input  signal-to-noise  ratio  for  the  v-th  law 
half-wave  detector  must  be  increased  by  the  factor 


F^(v)  =  I  (4D(v)  -  2)^^^  ,  (145) 

in  order  to  maintain  the  same  output  signal-to-noise  ratio.  This  factor  is 
tabulated  in  table  2  and  plotted  in  figure  9.  The  most  striking  feature  of 
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Table  2.  Additional  Input  Signal -to-Noi se  Ratio  Relative  to  Full-Wave 
Square-Law  System;  Low  Input  Signal-to-Noise  Ratio, 

Half-Wave  Rectifiers,  Selected  v,  Independent  Samples 
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these  results  is  the  large  degradation,  relative  to  the  full-wave  square-law 
system,  suffered  by  employing  a-  half-wave  rectifier  for  independent  samples. 
For  example,  a  linear  half-wave  rectifier  requires  an  additional  3.16  dB  input 
signal-to-  noise  ratio  in  order  to  realize  the  same  output  signal -to-noise 
ratio  as  a  full-wave  square-law  detector,  while  the  square-law  half-wave 
rectifier  requires  1.99  dB  additional  input  signal-to-noise  ratio. 

Amelioration  of  this  degradation  in  the  case  of  statistically  dependent 
samples  will  be  demonstrated  later  for  bandpass  spectra  with  various  values  of 
Q. 
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Figure  9.  Rdditional  Input  S i gn a  I -to-No i se 
Ratio  Relative  to  Square-Lauj  System;  Half-Wave 
Rectifiers,  Lou  Input  5 i gn a  1 -to-No i se  Ratio 
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FLAT  BANDPASS  SPECTRUM  AND  BOX  CAR  FILTER 

In  this  section,  the  input  spectrum  is  taken  to  be  flat  in  a  bandwidth 
W  about  center  frequency  +f^;  see  figure  3.  The  Q  of  the  spectrum  will  be 
varied  from  its  minimum  of  1/2  to  values  large  enough  that  the  limiting 
performance  achieved  by  a  narrowband  spectrum  is  virtually  achieved. 

The  sampled  version  of  the  detection  system,  depicted  in  figure  2,  will 
be  considered,  with  the  weights  {w(n)^  in  the  accumulator  set  equal  to  the 
same  value;  this  is  the  box  car  filter  of  (10).  Sampling  increment  a  is 
arbitrary,  as  is  the  power  law  v  of'  the  half-wave  and  full-wave  rectifiers  to 
be  considered.  When  a  is  taken  very  small,  the  system  performance  will 
approach  that  for  a  nonimpulsive  filter;  thus,  there  is  no  need  to 
additionally  evaluate  the  nonimpulsive  results  of  (95)  and  (103).  On  the 
other  hand,  when  a  is  taken  large,  the  effects  of  undersampling  will  become 
apparent. 

The  basis  of  comparison  for  the  u-th  law  rectifiers  considered  here  has 
been  delineated  in  an  earlier  section,  in  particular,  in  the  discussion 
surrounding  (59)  and  (60).  We  confine  the  numerical  results  for  the  decibel 
difference  in  (60)  to  low  input  signal -to-noi se  ratio,  and  specialize  the 
general  results  in  (91)  and  (102)  for  symmetric  full-wave  and  half-wave 
rectifiers,  accordingly.  The  detailed  computational  procedure  is  presented  in 
appendix  D,  for  both  the  symmetric  full-wave  rectifier  and  the  half-wave 
rectifier.  Sampling  increment  a  is  related  to  integration  time  or  observation 
time  T  according  to 
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T  =  MA  ,  (146) 

where  M  is  the  number  of  samples ' employed  in  box  car  filter  (10). 

Full-Wave  Rectifier;  Variable  Sampling  Increment 

In  this  subsection,  the  nonlinear  detector  is  taken  to  be  a  u-th  law 
full-wave  rectifier.  The  first  result  in  figure  10  gives  the  additional 
signal-to-noise  ratio  in  dB  required  by  the  full-wave  rectifier  relative  to 
the  standard  of  (57),  for  TW  =  50,  Q  =  1/2,  and  M  varied  from  50  up  to  500. 
Since  A  is  given  by  (146)  in  terms  of  integration  time  T  and  number  of  samples 
M,  there  follows,  for  the  sampling  frequency  f^  =  1/A,  the  ratio 

^  -  3_  _  M_ 

W  WA  ■  TW  • 

Thus,  for  example,  M  =  100  here  corresponds  to  a  ratio  of  sampling  frequency 
f^  to  bandwidth  W  of  100/50  =  2;  since  Q  =  1/2,  this  value  of  M  corresponds 
to  sampling  of  an  input  process  with  a  lowpass  spectrum,  at  twice  the  highest 
frequency,  meaning  independent  samples.  These  frequencies  are  indicated  in 
the  figure.  The  dB  numbers  entered  at  this  value  of  M  on  figure  10  for  v  =  1, 
2,  3,  are  in  fact  exactly  those  already  listed  earlier  in  table  1  for 
independent  samples. 

Increasing  M  to  200  leads  to  a  sampling  frequency  4  times  f^(=  W) ,  and 
to  dB  values  virtually  equal  to  the  saturation  values  entered  at  M  =  500. 

Thus,  it  is  possible  to  lower  the  required  input  signal-to-noise  ratio  by 
.133,  .018,  .077  dB  for  u  =  1,  2,  3,  respectively,  by  employing  larger  M 
values  than  lead  to  independent  samples.  The  slightly  negative  dB  values  for 
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Figure  10.  Flat  Spectrum,  Box  Car  Filter, 
Full-Wave  Rectifier,  TW  =  50,  0  =  1/2 
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Figure  11.  Flat  Spectrum,  Box  Car  Filter, 
Full-Wave  Rectifier,  TW  =  50,  Q  =  5/6 
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large  M,  achieved  by  the  full-wave  rectifier  with  v  =  2,  are  consistent  with 
the  discussion  surrounding  (58),' and  merely  reflect  the  fact  that  the  standard 
is  not  the  absolute  optimum  output  signal-to-noise  ratio. 

On  the  other  hand,  decreasing  the  sampling  frequency  below  2f^,  that 
is,  decreasing  M  below  100  for  this  example  of  TW  =  50,  leads  to  a  significant 
degradation  of  performance,  as  witnessed  by  the  rapid  rise  of  the  dB  curves  in 
figure  10,  to  the  left  of  M  =  100. 

In  figure  11,  the  only  change  is  to  increase  Q  to  5/6;  this  is  a  flat 
bandpass  spectrum.  Since,  in  general,  using  (12)  and  (146),  the  ratio 


TW 


(148) 


f^  AW(1  +  20) 


then  the  value  of  M  =  133  in  figure  11  corresponds  to  f^  =  2f^  =  (8/3)W. 

Again,  this  is  seen  to  correspond  to  a  prominent  knee  of  the  performance 
degradation  curve.  Increasing  M  to  250  realizes  essentially  the  same  input  dB 
values  as  for  M  =  500;  however,  for  u  =  2,  M  need  only  be  increased  to  150  in 
order  to  essentially  attain  the  large-M  asymptote.  It  is  worthwhile  noting 
for  future  reference  that  the  frequency  ratios  in  (147)  and  (148)  are  directly 
proportional  to  M,  the  number  of  samples  employed  in  the  box  car  filter.  The 
seemingly  anomalous  behavior  of  figure  11  near  M  =  75  is  more  pronounced  in 
the  next  figure  and  will  be  explained  there. 

In  figure  12,  the  only  change  is  to  increase  Q  to  2,  while  keeping 
TW  =  50.  Reference  to  (147)  and  (148)  reveals  that  M  =  250  now  corresponds  to 
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Figure  12.  Flat  Spectrum,  Box  Car  Filter, 
Full-Wave  Rectifier,  TW  =  50,  0=2 
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=  2f^  =  5W;  there  is  a  slight  improvement  by  increasing  M  to  300. 

However,  the  most  striking  feature  of  this  figure  is  the  return  of  the 
required  input  signal-to-noise  ratio  dB  values  to  very  low  levels  for  M  in  the 
range  (125,150).  To  explain  this  behavior,  we  first  note  from  (147)  and 

(148)  that  f^  =  f^  =  2.5W  at  M  =  125,  while  f^  =  1.2f^  =  3W  at  M  =  150;  thus, 
the  sampling  rate  is  still  larger  than  twice  the  bandwidth  W  of  the  input 
process  (and  remains  so  until  M  decreases  to  100).  However,  since  f^  is 
smaller  than  2f^,  we  have  an  example  of  undersampling,  without  significant 
loss  in  performance. 

The  mathematical  explanation  of  what  is  happening  requires  a  close 
examination  of  the  denominator  of  the  system  output  signal-to-noise  ratio, 

Yjp,  in  (91).  In  appendix  E,  some  useful  frequency  domain  representations 
of  the  denominator  of  (91)  are  derived.  In  particular,  a  combination  of 
(E-8)  and  (E-12)  yields  the  approximate  result 


(149) 


where 


for  all  f 


(150) 


n 


(k) 

IS  the  aliased  (periodic)  version  of  P  (f),  which,  in  turn,  is  the  Fourier 

k 

transform  of  p  (T);  see  (47).  The  approximation  (149)  is  valid  when  the  TW 


product  is  large;  see  (E-8)  for  details. 
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The  leading  term  in  the  denominator  of  signal-to-noise  ratio  in 
(91)  is  equal  to  U(2)  times  the  result  (149).  The  latter  quantity  is  the 
average  value  under  the  square  of  the  aliased  input  spectrum.  In  order  to 
realize  large  values  of  y^p,  then  (149)  should  be  small.  Plots*  of 
P^^^(f)  in  figures  13  and  14  for  Q  =  2,  with  M  =  125  and  M  =  150, 
respectively,  reveal  that  none  of  the  aliased  lobes  overlap,  despite  the 
undersampling;  the  cross-hatched  lobes  represent  the  input  (unaliased) 
spectrum.  Thus,  the  integral  of  (149)  remains  constant  (and  small)  for  M  in 
the  range  (125,150) . 

However,  if  M  is  decreased  below  125,  the  aliasing  lobes  in  figure  13 
bordering  f  =  0,  for  example,  begin  to  overlap,  thereby  doubling  the  value  of 
‘P^^^(f)  in  this  region.  Similarly,  if  M  is  increased  above  150,  the  lobes 
in  figure  14  bordering  f/f^  =  +.75  begin  to  overlap.  In  both  these  cases, 
the  value  of  (149)  would  increase  significantly,  thereby  decreasing  the  system 
output  signal-to-noise  ratio  y^p  in  (91),  or  alternatively  increasing  the 
required  input  signal-to-noise  ratio.  In  order  to  avoid  this  degradation, 
overlapped  aliasing  lobes  must  be  avoided,  either  by  looking  for  "clean" 
regions  for  high  Q,  as  in  figures  13  and  14,  or  by  sampling  at  rates  greater 
than  twice  the  highest  frequency  of  the  input  spectrum. 


*A  slight  tilt  has  been  added  to  the  spectral  shape  in  order  to  more  easily 
identify  which  lobes  result  from  positive  and  negative  frequency  components. 
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The  results  in  figure  12,  as  well  as  all  the  other  figures,  have  not  been 
limited  to  considering  just  the  k  =  2  term  in  the  denominator  of  the  system 
output  signal-to-noi se  ratio  y,  nor  have  they  utilized  the  approximations 
(149)  or  (E-8).  Instead,  they  have  been  obtained  by  utilizing  K  =  1000  terms 
in  the  series,  and  by  actually  conducting  the  exact  summations  in  the 
denominators  of  (91)  and  (102).  The  discussion  above,  relative  to  figures  13 
and  14,  was  presented  in  order  to  give  a  simple  physical  explanation  of  what 
is  happening,  and  thereby  furnish  guidance  to  further  cases  of  interest. 

In  figure  15,  Q  is  further  increased  to  3,  again  keeping  TW  at  50. 

Reference  to  (147)  and  (148)  indicates  that  M  =  350  corresponds  to 

f^  =  2f^  =  7W;  slightly  improved  performance  can  be  achieved  by  increasing 

M  to  400.  Substantially  the  same  performance  level  can  be  achieved,  however, 

for  M  in  the  range  (175,250),  corresponding,  respectively,  to  undersampling 

with  f^  in  the  range  3.5W  to  5W.  Furthermore,  there  is  an  additional 

possibility  for  M  in  the  range  (117,125),  with  sampling  frequencies  in  the 

range  f^  =  (2.34W,  2.5W);  it  may  be  verified  by  use  of  (147)  and  (148)  that 

1  ^ 

the  aliasing  lobes  of  P'  '(f)  do  not  overlap  for  this  range  of  sampling 
frequencies  (as  shown  in  figures  13  and  14).  However,  only  the  v  =  2 
full-wave  rectifier  achieves  the  level  attainable  at  large  M,  whereas  the 
u  =  1  and  u  =  3  rectifiers  suffer  additional  degradation.  This  is  due  to  the 
fact  that  the  denominator  series  in  (91)  terminates  for  v  =  2,  but  does  not 
terminate  for  v  =  1  or  3;  thus,  the  higher-order  terms  in  (91)  cause 
additional  unavoidable  degradation  through  leakage  of  higher-order  spectral 
terms,  as  indicated  in  appendix  E. 
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Number  of  Samples,  M 

Figure  15.  Flat  Spectrum,  Box  Car  Filter, 
Full-Wave  Rectifier,  TW  =  50,  0=3 
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Figure  16.  Variation  with  Q,  Full-Wave 
Reotifier,  TW  =  50,  M  =  1000 
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Full-Wave  Rectifier;  Variable  Q 

The  curves  in  the  previous  subsection  all  saturated  for  large  numbers  of 
samples,  M;  that  is,  the  required  input  signal-to-noise  ratio  approached  that 
for  a  continuous  box  car  filter.  Here,  we  will  take  M  =  1000,  meaning  that  we 
are  effectively  considering  only  the  continuous  box  car  filter,  and  look  now 
at  the  required  system  input  signal-to-noise  ratio  as  a  function  of  the  Q  of 
the  spectrum,  but  not  limited  to  the  four  discrete  values  of  1/2,  5/6,  2,  and 
3  earlier.  The  first  result  in  figure  16  indicates  a  rapid  drop  in  required 
input  signal-to-noise  ratio  as  Q  increases  above  .5.  The  region  (.5,. 6)  for  Q 
is  blown  up  in  figure  17  to  better  illustrate  the  decay.  The  main  conclusion 
from  these  two  figures  is  that  once  Q-.5  is  larger  than  approximately  1,  the 
required  input  signal-to-noise  ratio  is  essentially  independent  of  Q. 

To  illustrate  the  relative  independence  of  the  results  on  the  specific  TW 
product  (once  it  is  large),  figure  18  was  computed  and  plotted  for  TW  =  100; 
this  result  can  be  compared  with  figure  17  which  utilized  TW  =  50. 

Finally,  in  a  similar  comparison,  figure  11  for  TW  =  50,  Q  =  5/6  was 
rerun  for  TW  =  100,  Q  =  5/6  in  figure  19.  Except  for  the  doubling  of  the 
abscissa  values  (M  =  1000  vs  500  earlier),  figures  19  and  11  are  very 
similar.  The  fact  that  M  must  be  doubled  if  TW  is  doubled  is  reasonable, 
whether  that  doubling  comes  about  from  increased  T  (more  observation  time)  or 
increased  W  (higher  bandwidth,  meaning  faster  sampling). 
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Qual i ty  F actor ,  Q 

Figure  17.  Small  Q  Variation,  Full-Wave 
Rectifier,  TW  =  50,  M  =  1000 


Figure  18.  Small  Q  Variation,  Full-Wave 
Rectifier,  TW  =  100,  M  =  1000 
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Figure  19.  Flat  Spectrum,  Box  Car  Filter, 
Full-Wave  Rectifier,  TW  =  100,  Q  =  5/6 


Number  of  Samples,  M 

Figure  20.  Flat  Spectrum,  Box  Car  Filter, 
Half-Wave  Rectifier,  TW  =  50,  Q  =  1/2 
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Half-Wave  Rectifier;  Variable  Sampling  Increment 

We  will  repeat  the  series  of  plots  given  earlier  in  figures  10,  11,  12, 
15,  but  now  for  w-th  law  half-wave  rectifiers.  The  first  result  in  figure  20 
gives  the  required  input  signal-to-noise  ratio  (relative  to  the  standard)  for 
0  =  1/2,  a  low-pass  spectrum.  The  levels  of  required  signal-to-noise  ratios 
are  much  greater  here  for  the  half-wave  rectifiers  than  for  the  corresponding 
full-wave  rectifiers  in  figure  10.  M  =  100  corresponds  to  sampling  frequency 
f^  =  2f^  =  2W;  increasing  M  to  200  essentially  reaches  the  saturation 
value  attained  for  large  M  (continuous  filtering).  The  dB  numbers  entered  at 
M  =  100  agree  with  those  given  earlier  in  table  1,  since  this  corresponds  to 
independent  samples.  (When  TW  is  doubled  to  the  value  100,  virtually  the  same 
plot  results  when  the  abscissa  scale,  M,  is  also  doubled.) 

When  0  is  increased  to  5/6  in  figure  21,  a  marked  improvement  in 
performance  occurs.  The  value  M  =  133  corresponds  to  f^  =  2f^  =  (8/3)W; 
increasing  M  to  200  achieves  saturation  values,  at  which  point  f^  =  3f^  =  4W. 
However,  the  performance  relative  to  the  corresponding  full-wave  rectifier 
result  in  figure  11  is  still  poorer. 

For  0  =  2  in  figure  22,  although  there  is  a  marked  dip  in  the  required 
input  signal-to-noise  ratio  for  M  in  the  range  (130,150),  the  levels  achieved 
are  not  as  low  as  those  possible  for  larger  M.  For  example,  an  additional 
.24  dB  is  required  at  M  =  150  for  v  =  1  than  at  M  =  400  for  u  =  2.  At 
M  =  125,  the  dB  values  in  figure  22  are  very  large;  this  is  due  to  an  aliased 
component,  a  la  figures  13  and  14,  abutting  frequency  f  =  0.  This  result  is 
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Figure  21.  Flat  Spectrum,  Box  Car  Filter, 
Half-Wave  Rectifier,  TW  =  50,  Q  =  5/6 
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Figure  22.  Flat  Spectrum,  Box  Car  Filter, 
Half-Wave  Rectifier,  TW  =  50,  0=2 
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distinctly  different  from  the  corresponding  one  for  a  full-wave  rectifier  in 
figure  12. 

The  half-wave  rectifier  result  for  Q  =  3  in  figure  23  is  somewhat  similar 
to  the  corresponding  full-wave  rectifier  result  in  figure  15.  However,  there 
is  much  more  degradation  here  for  M  near  120,  and  moderate  degradation  for  M 
in  the  range  (180,250).  Essentially  the  same  performance  for  u  =  2  is 
achieved  at  M  =  350,  f^  =  2f^  as  at  M  =  500,  f^  =  (20/7)f^,  but  losses 
are  encountered  for  sampling  rates  between  these  values. 

Half-Wave  Rectifier:  Variable  0 

To  determine  the  explicit  dependence  of  performance  on  the  Q  of  the  input 
spectrum,  we  eliminate  the  dependence  on  M,  by  setting  M  equal  to  a  large 
value,  namely  1000,  thereby  essentially  realizing  a  continuous  box  car 
filter.  Then  we  vary  Q  over  the  range  (.5,3)  and  plot  the  results  for  TW  =  50 
in  figure  24.  Saturation  is  achieved  for  Q  >  1.5  at  acceptably  low  values; 
however  the  performance  degrades  considerably  for  small  Q. 

Comparison  of  figure  24  with  the  corresponding  result  in  figure  16  for  a 
full-wave  rectifier  (both  with  TW  =  50  and  M  =  1000)  reveals  that  by  the  time 
0  reaches  1.5,  the  half-wave  and  full-wave  rectifiers  realize  virtually  the 
same  performance  level,  regardless  of  the  value  of  v.  For  example,  the  dB 
numbers  listed  above  Q  =  3  are  almost  identical.  This  confirms  the  results 
anticipated  in  the  analysis  presented  earlier  in  (78)-(80)  et  seq. 
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Figure  23.  Flat  Spectrum,  Box  Car  Filter, 
Half-Wave  Rectifier,  TW  =  50,  Q  =  3 


Figure  24.  Variation  with  Q,  Half-Wave 
Rectifier,  TW  =  50,  M  =  1000 
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The  region  (.5,. 6)  for  Q  is  blown  up  in  figure  25;  it  indicates 
approximately  that  if 


a  plateau  in  performance  is  essentially  realized,  while  a  more  gradual 
improvement  takes  place  for  larger  Q.  This  same  bound  on  Q  was  previously 
encountered  in  (111)  when  we  discovered  the  condition  under  which  the  linear 
half-wave  rectifier  and  the  symmetric  linear  full-wave  rectifier  have  similar 
performance.  The  analysis  in  (109)-(117)  is  again  directly  relevant  here  in 
regards  to  overlapping  aliased  spectral  lobes. 

If  TW  is  increased  to  100,  figure  26  illustrates  very  similar  behavior, 
except  that  the  sharp  transition  region  near  Q  =  1/2  is  approximately  half  as 
wide;  this  is  consistent  with  (151). 

In  a  similar  vein,  for  a  fixed  Q  of  5/6,  a  comparison  of  results  for 
TW  =  100  in  figure  27  can  be  made  with  the  corresponding  earlier  result  for 
TW  =  50  in  figure  21.  Except  for  the  doubling  of  the  M  scale  in  figure  27, 
the  two  sets  of  curves  are  virtually  identical. 
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Figure  25.  Small  Q  Variation,  Half-Wave 
Rectifier,  TW  =  50,  M  =  1000 
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Figure  26.  Small  Q  Variation,  Half-Wave 
Rectifier,  TW  =  100,  M  =  1000 
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Figure  27.  Flat  Spectrum,  Box  Car  Filter, 
Half-Wave  Rectifier,  TW  =  100,  Q  =  5/6 
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OTHER  SPECTRA  AND  FILTERS 


Gaussian  Spectrum 


All  the  previous  results  pertained  to  the  flat  bandpass  spectrum  of 
figure  3.  In  order  to  ascertain  how  important  the  details  of  the  input 
spectrum  are,  we  consider  in  this  subsection  the  input  normalized  correlation 


p{T)  =  cos(2irf^T)  exp 


(-  2  • 


(152) 


with  corresponding  spectrum  given  by  (47)  as 


2^^^W 


exp 


-2ir 


if^ 

W 


+  exp  - 


(153) 


This  is  a  pair  of  Gaussian  lobes  centered  at  +f 


It  should  be  observed  that 


(154) 


independent  of  center  frequency  f^.  Thus  W  in  (152)  and  (153)  is  the 
effective  or  statistical  bandwidth  of  the  positive  lobe  of  input  spectrum 
P^^\f),  as  if  the  negative  lobe  were  absent.  Also,  the  relative  value  of 
the  positive  lobe  at  frequencies  f  =  f^.  ±  W/2  is  exp(-ir/2)  =  .208  =  -6.82  dB. 
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More  generally,  the  effective  bandwidth  of  the  complete  spectrum 


P^^^f)  is 


ffdf  P^^^(f)]^  ^  _ 1 

r  r 


Jdf  P^^^  (f)  Jdf  P^^^  (f) 


1  +  exp(-4irQ  ) 


2W 


2 


(155) 


Unfortunately,  this  depends  on  Q  =  fj./W.  If  we  were  to  use  this  latter 
result  in  standard  output  signal-to-noise  ratio  in  (51),  we  would  have 

a 


(156) 


1  +  exp(-4ir0  ) 


and  (56)  would  yield 


(157) 


The  dependence  on  Q  is  undesirable,  although  it  is  a  weak  dependence;  for 
example,  for  Q  =  1/2,  the  exponential  in  (157)  is  .043.  Hence  we  drop  this 
dependence,  and  use  the  usual  basis  (57)  again. 

If  we  use  (91)  for  or  (102)  for  y  ,  with  the  S+N  versus  N 

I  r  in 

hypotheses  of  (56),  and  a  small  input  signal-to-noise  ratio,  we  have 


^I(F,H)  2 


(158) 


where 


(159) 


k 


n 
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the  sum  being  over  k  =  2,  4,  6,  ...  for  full-wave  rectifiers,  and  over 
k  =  1,  3,  5,  ...  for  half-wave  rectifiers.  If  we  now  equate  the  output 
signal-to-noise  ratios,  that  is,  set 
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(160) 


then  there  follows  factor 


$/N  2 


1/2 


(S/N)^  V 


(161) 


This  is  the  input  signal-to-noise  ratio  required,  relative  to  basis  (57). 
Recall  that  W  is  now  the  effective  BW  of  the  positive  frequency  lobe;  see 
(154). 

Plots  of  the  relative  increase  in  input  signal-to-noise  ratio,  as  given 
by  (161)  and  (159),  are  presented  in  figures  28  and  29  for  the  full-wave  and 
half-wave  rectifiers,  respectively,  for  TW  =  50  and  Q  =  5/6.  These  should  be 
compared  with  the  corresponding  flat  bandpass  spectrum  results  in  figures  11 
and  21,  respectively.  Except  for  a  general  smoothing  in  figures  28  and  29, 
due  to  the  smoother  Gaussian  spectrum  (153),  the  results  are  very  similar. 
Saturation  is  essentially  reached  at  M  =  200,  which  corresponds  to  sampling 
frequency  f^  =  3f^  =  4W. 
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Figure  28.  Gaussian  Spectrum,  Box  Car  Filter, 
Full-Wave  Rectifier,  TW  =  50,  Q  =  5/6 


Figure  23.  Gaussian  Spectrum,  Box  Car  Filter, 
Half-Wave  Rectifier,  TW  =  50,  Q  =  5/6 
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-6  dB/Octave  Spectrum 

2 

In  this  section,  the  input  spectrum  is  taken  to  have  a  1/f  shape  in 
the  band  There  is  a  difference  in  this  example  between  the 

bandwidth  W  =  f2  -  f-j  and  the  effective  bandwidth,  W^,  which  is 
developed  in  appendix  F.  The  pertinent  equations  for  this  case  (as  well  as 
the  spectral  shapes  f^  for  n  =  -2, -1,0, 1,2)  are  also  presented  in  appendix 
F.  The  results  for  TW  =  50  and  Q  =  5/6  are  given  in  figures  30  and  31  for 
full-wave  and  half-wave  rectifiers,  respectively.  They  are  very  similar  to 
the  corresponding  earlier  results  with  the  same  parameter  values. 

RC  Filter 


All  the  previous  examples  have  employed  a  box  car  impulse  response,  as 
given  by  (10).  We  now  replace  this  assumption  by  one  in  which  the  impulsive 
filter  weights  are  given  by  samples  of  an  RC  filter  response,  with  effective 
duration  T  =  2  RC;  see  (52)  or  (72).  That  is. 


w(n)  =  ^  exp 


for  n  >  0 


(162) 


The  effective  number  of  samples  is,  from  (67), 


M  = 


^  w(n) 


n 


1  +  exp 


2^/, 


12  R^C 


f  RC  ^  , 

2.2  J  —  >  1 


(163) 
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Figure  30.  -6  dB/Octave  Spectrum,  Box  Car 
Filter,  Full -Wave  Rectifier,  TW  =  50,  Q  =  5/6 


Number  of  Samples,  M 

Figure  31.  -6  dB/Octave  Spectrum,  Box  Car 
Filter,  Half-Wave  Rectifier,  TW  =‘50,  Q  =  5/6 
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For  example, 


.  RC 
for  —  = 
A 


10  .  M  s  20 


120o)  • 


Accordingly  we  set 


pop  T 

M  =  j  (for  M  »  1)  .  (164) 

This  is  consistent  with  (146)  employed  for  the  box  car  filter.  The  normalized 
autocorrelation  is,  via  (68), 


B(n)  =  exp 


(165) 


A  plot  of  the  relative  input  signal-to-noise  ratio  is  given  in  figure  32 
for  a  half-wave  rectifier  with  TW  =  50,  Q  =  5/6,  and  a  flat  bandpass  input 
spectrum.  It  indicates  the  same  general  behavior  as  corresponding  earlier 
results  for  different  input  spectra  and/or  filters;  the  precise  numerical 
values  are  a  little  different,  and  are  a  reflection  of  the  particular  filter 
employed  here,  namely  (162). 
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Figure  32.  Flat  Spectrum,  RC  Filter, 
Half-Wave  Rectifier,  TW  =  50,  Q  =  5/6 
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A  PITFALL  OF  THE  LONG  AVERAGING  TIME  ASSUMPTION 


The  general  factor  by  which  the  input  signal-to-noise  ratio  for  a 
half-wave  u-th  law  rectifier  must  be  increased,  relative  to  the  standard 
full-wave  square-law  rectifier,  is  given  in  (D-5).  Here,  Sum^^  is  given  in 
(0-1),  in  terms  of  the  normalized  autocorrelation  (68)  of  the  filter  weights. 
Under  the  long  averaging  time  assumption,  as  discussed  earlier  with  regard  to 
(74)  et  seq.,  the  quantity  B(n)  is  replaced  by  its  origin  value  of  1  for  all 
n.  Then  the  quantity  S(k)  in  (0-3)  is  replaced  by 


oo 


PQ(nA)  for  k  >  1  , 


(166) 


n=l 


and  factor  (0-5)  is  replaced  by 


(167) 


where 


(168) 


k=l 


is  an  obvious  modification  of  (0-2). 


A  similar  approach  for  full-wave  w-th  law  rectifiers  yields 
signal-to-noise  ratio  factor 


(169) 


where 
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oo 

Sp  =  D(v)  -1+2  ^  U(l<)  S(k)  .  (170) 

k=2 

even 

It  will  be  seen  from  (168)  and  (170)  that  a  fundamental  calculation  required 
is  that  described  in  (166),  whether  utilizing  the  full-wave  or  the  half-wave 
rectifier.  In  particular,  for  the  flat  bandpass  input  spectrum  given  by 
figure  3,  with  normalized  correlation  (13),  the  summation  in  (166)  is  very 
slowly  convergent  for  low  order  k.  This  problem  is  treated  in  appendix  G, 
through  a  judicious  combination  of  Poisson's  formula  and  numerical  calculation. 

Results  for  the  long  averaging  time  assumption  are  superposed  as  dotted 
lines  in  figures  33  and  34  over  the  earlier  results  from  figures  24  and  25, 
respectively,  for  half-wave  rectifiers,  TW  =  50,  M  ==  1000.  As  Q  approaches 
.5,  there  is  a  marked  difference  between  the  two  results.  For  example,  for 
V  =  2,  the  approximation  yields  .34  dB  whereas  the  exact  results  is  1.93  dB. 

And  for  V  =  1,  the  comparison  is  .17  dB  versus  exact  value  3.12  dB,  almost  a  3 
dB  discrepancy.  This  serves  to  point  out  the  pitfall  of  employing  the  long 
averaging  time  assumption  when  inappropriate.  Regardless  of  the  size  of  the 
TW  product,  there  will  always  be  a  narrow  range  near  Q  =  1/2  where  the  exact 
signal-to-noise  ratio  drops  sharply  from  the  values  listed  in  tables  1  and  2 
to  their  eventual  asymptotic  values  for  large  Q.  This  transition  is  ignored 
by  the  long  averaging  time  assumption  and  can  be  overly  optimistic  in  its 
performance  prediction. 
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Figure  33.  Long  Rveraging  Time  Assumption, 
Half-Nave  Rectifier,  TN  =  50,  M  =  1000 


Quality  Factor,  0 


Figure  34.  Long  Rveraging  Time  Assumption, 
Half-Nave  Rectifier,  TN  =  50,  M  =  1000,  Small  Q 
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SUMMARY 

The  losses  incurred  by  using  various  full-wave  and  half-wave  v-th  law 
rectifiers  with  different  filter  characteristics,  sampling  rates,  and  input 
spectra  have  been  evaluated  and  plotted  for  a  number  of  cases.  These  results 
allow  for  a  ready  comparison  of  many  alternative  choices  and  give  quantitative 
bases  for  a  selection  procedure. 


The  possibility  of  deliberately  undersampling  an  input  process  with  a 
high-Q  spectrum,  with  insignificant  loss  of  performance,  has  been  analyzed  and 
numerically  investigated.  The  exact  choice  of  sampling  frequency  is  crucial, 
but  can  be  easily  calculated  once  the  lower  and  upper  frequencies  of  the 
band-limited  input  spectrum  are  specified.  Even  though  the  system 
nonlinearity  g  [x]  creates  harmonics  and  intermodulation  products,  and  the 
undersampling  process  creates  aliased  spectral  replicas,  it  is  still  possible, 
through  proper  choice  of  the  sampling  rate,  to  control  all  these  undesirable 
by-products  and  achieve  a  near-minimum  input  signal-to-noise  ratio  for 
specified  output  deflection. 


The  danger  of  using  the  long  averaging  time  assumption  when  inappropriate 
has  been  demonstrated  via  numerical  example.  All  the  results  presented  in 
this  report  have  not  employed  this  assumption,  but  have  utilized  the  precise 
filter  characteristics  and  finite  time  extent.  Programs  for  these  evaluations 
are  presented  in  tables  H-1  through  H-4  in  appendix  H.  Changes  required  to 
handle  the  Gaussian  spectrum  or  the  RC  filter,  instead  of  the  flat  bandpass 
spectrum  or  box  car  filter,  respectively,  are  also  presented  there. 
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The  losses  incurred  by  a  half-wave  rectifier  for  low-0  inputs  are  very 
significant  and  should  be  avoided,  either  by  utilizing  a  full-wave  rectifier 
or  by  filtering  out  the  low  frequency  components  of  the  input,  prior  to 
sampling  and  nonlinear  distortion.  Generally  speaking,  the  square-law 
detector,  v  =  2,  offers  the  best  performance. 

Direct  comparisons  of  these  results  with  those  in  [1]  are  not  possible, 
because  there  are  no  common  examples.  Also,  most  of  the  results  in  [1]  are 
for  very  narrowband  spectra  and  employ  a  long  averaging  time  assumption,  in 
addition  to  using  very  few  terms  in  the  series  expansion  for  the  output 
variance. 
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APPENDIX  A.  DERIVATION  OF  RECURRENCE 


The  quantity  of  interest  is  given  by  (85): 

CO 

L(k)  =  J'  dw  w’’  0(w)  He|^(w)  for  ail  k  >  0  . 
o 


Fi  rst  of  an  ,  using  (21 ) , 
L(0)  =  J  dw  w’' 


QO  00 

0(w)  =  (2ir)"^^^  dw  w’'  exp(-w^/2)  = 
o 


J-1  -1 


w  w£L 

=  (2ir)"^^^  J  dx  (2x)  ^  exp(-x)  =2^  ir  ^ 
0 

In  a  similar  fashion,  there  follows 


r("f) 


(A-l) 


(A-2) 


«-l  1 


L(l) 


=  I  dw  w^-"^  0(w)  =  2  2  2  p^i  ^ 


(A-3) 


We  now  employ  the  recurrence  [3,  (22.7.14)] 


He|^(w)  =  w  Hej^_.|(w)  -  (k-1)  He|^_2(w)  for  k  >  2 


(A-4) 


Substitution  in  (A-l)  yields 


L(k)  =  ^dw  w’^"'’^  0(w)  Hej^_.|(w)  -  (k-1)  L(k-2)  for  k  >  2 


We  now  integrate  by  parts,  using 


(A-5) 


.v+1 


u  =  w  ,  dv  =  dw  ia(w)  He|^_^(w) 


du  =  dw  (u+1)  w  ,  V  =  -0(w)  Hej^_2(w)  , 


(A-6) 
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the  last  relation  from  [3,  (22.11.8)].  The  integral  in  (A-5)  then  becomes, 
since  w  >  0, 


0(w)  He|^_2(w)J^  +  (v+1) 

=  (v+1)  L(k-2)  . 


dw  w’’  0(w)  He|^_2(w)  = 

0 


(A-7) 


Use  of  this  result  in  (A-5)  immediately  yields 

L(k)  =  (v+2-k)  L(k-2)  for  all  k  >  2  .  (A-8) 

The  starting  values  for  this  recurrence  have  already  been  furnished  in 
(A-2)  and  (A-3). 
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APPENDIX  B.  DERIVATIONS  OF  (110)  AND  (112) 


We  employ  (73)  and  (48)  to  develop  the  left-hand  side  of  (110)  according 
to 


dr  a(r)  p(T)  =  dr 


du  h(u)  h(u-'c) 


du  h‘^(u) 


df  exp(i2irfT)  P^^\f)  = 


Ijdu  h^(u)J 


-1 


df  P^^^f)  Idu  h(u) 


dr  h(u-T)  exp(i2irfr)  = 


=  Idf  P^^^f)  lH(f)l^ 

Jdf  |H(f)l  ^ 


(B-1) 


where  we  let  x  =  u-t  in  the  innermost  r  integral  in  the  second  line,  and  used 
Parseval's  theorem  for  the  denominator  term. 


For  the  impulsive  filter  case,  we  utilize  (68)  and  (48)  to  express  the 
left-hand  side  of  (112)  as 


^  B(n)  p(na) 
n 


2 


^  w(m)  w(m-n) 

JT] _ 

5-  w^(m) 
m 


J 


^df  exp(i2irfan)  P^^^f)  = 


^  w^(m) 
.  m 


-1 


w(m-n)  exp(i2irfAn)  = 


r  df  P^^^f)  |H(f)l  ^ 

^  w^(m) 


(B-2) 
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where  we  let  k  =  m-n  in  the  innermost  sum  on  n  in  the  second  line,  and  used 

H(f)  =  ^  w(n)  exp(-i2irfan)  ,  (B-3) 

n 

which  follows  from  (5)  and  (8),  Since  H(f)  obviously  has  period  1/A  in  f,  an 
integral  over  any  period  of  the  power  transfer  function  becomes 


r  df  lH(f)l  ^  = 


1 

A 


2  w^Cm)  =  \  b(0) 
m 


(B-4) 


using  orthogonality  of  the  exponentials  for  m  ;>«  n.  Thus  (B-2)  translates  into 


^  I3(n)  p(nA) 
n 


fdf  P^^^fl  iH(f)l  ^ 

A  r  df  |H(f)|^ 


Equation  (B-4)  is  a  special  case  of  the  useful  result  that 


Jdf  lH(f)p  exp(i2irfAn)  =  \^  w(m) 

'/a  m 


w(m-n)  =  -  b(n) 


(B-5) 


(B-6) 


The  method  of  derivation  is  similar  to  (B-4),  and  (62)  has  been  employed  for 
the  final  identification. 
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APPENDIX  C.  LIKELIHOOD  RATIO  PROCESSOR  FOR  INDEPENDENT  SAMPLES 


For  a  zero  mean  Gaussian  input  process  with  variance  ,  the  joint 
probability  density  function  of  M  independent  samples  x. ,  ....  x„  is 

I  M 


Mr  /  2 

X. 


. a!' 


(c-1) 


The  likelihood  ratio  for  hypothesis  H-|  versus  Hq  is 


Pj^^l . ^M^  ^  l°S 

Po^^l . ^M^  VI 


exp 


V'^O  '’l  /  m=l 


m 


(C-2) 


The  sufficient  statistic  of  the  samples  is  obviously 


m  ’ 


m=l 


(C-3) 


which  is  interpreted  as  an  equi-weighted  sum  of  the  squared  input  samples. 
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APPENDIX  D.  COMPUTATIONAL  PROCEDURES 


The  denominator  of  output  signal-to-noise  ratio  in  (102)  for 

1 H 

half-wave  rectifiers  is 


OO 

2 

U(k)  2  B(n)  p^(nA)  = 

^  U(k) 

1  +  2 

cO  -j 

2.  B(n)  Po(nA)t  (0-1) 

k=l 

n 

k=l 

L. 

n=l 

=  2D(v)  -  1  +  2  ^  U(k)  S(k)  ,  (D-2) 

k=l 


where  we  used  (141)  and  defined  the  inner  sum  on  n  in  (D-1)  as 

oO 

S(k)  =  ^  B(n)  pQ(nA)  for  k  >  1  .  (0-3) 

n=l 

The  decay  of  S(k)  with  k  is  fast  enough  that  the  infinite  sum  on  k  in  (D-2) 
can  be  terminated  and  yet  realize  very  accurate  results  for  the  quantity 
Sum^. 


For  the  S-t-N  versus  N  input  hypotheses  of  (131)  and  low  input 
signal-to-noise  ratio,  the  output  signal-to-noise  ratio  (102)  becomes 
approximately 


^IH 


Sum^^  4 


for 


(0-4) 


If  this  is  equated  to  the  standard  value  in  (57)  with  input  signal-to-noise 
ratio  (S/N)^,  the  amount  by  which  the  input  signal-to-noise  ratio  for  the 
half-wave  rectifier  must  be  increased  is  the  factor 
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S/N  _  2 
(S/N)^  V 


(D-5) 


This  quantity  is  plotted  in  dB  for  a  variety  of  values  of  w,  sampling  rates, 
observation  times,  and  input  spectra,  as  discussed  in  (59)  and  (60).  The  way 
in  which  sampling  increment  A  is  related  to  the  effective  duration  T  of  the 
filter  is  according  to 


T  =  MA  ,  (D-6) 

where  M  is  the  effective  number  of  samples  of  the  filter.  Thus  the  ratio  of 
sampling  increment  to  inverse  bandwidth  is 


A 

1/W 


(D-7) 


in  terms  of  the  time-bandwidth  product  TW  and  number  of  samples  M. 


For  the  full-wave  rectifier,  the  only  change  is  to  limit  the  summation  in 
(D-1)  to  even  k,  with  the  result 


Sumr 


oo  <o 

=  2  U(k)  1  +  2  ^  (3(n)  PQ(nA) 


k=2 

even 


n=l 


J 


=  D(v)  -  1  +  2  ^  U(k)  S(k)  , 

k=2 

even 


(D-8) 


where  we  used  (130)  and  (D-3).  The  factor  corresponding  to  (D-5)  is  now 


S/N  ^  2 
(S/N)^  V 


(D-9) 


for  low  input  signal-to-noise  ratio.  This  quantity  has  been  plotted  in  dB  for 
various  v,  TW,  M,  and  input  spectra. 
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APPENDIX  E.  EQUIVALENT  FREQUENCY  DQMAIN  REPRESENTATIONS 


In  this  appendix,  we  derive  some  useful  equivalent  representations  of  the 
denominator  terms  in  the  output  signal-to-noise  ratios  of  (91)  and  (102).  We 
begin  by  considering  a  general  function  g(T)  with  Fourier  transform 


G(f)  = 


I 


dF exp(-i2irfT)  g(x)  . 


(E-1) 


Then  the  summation  of  interest  is 


^  B(n)  g(na)  =  ^  J  lH(f))  ^  exp(-i2irf An)  = 

n  "  .  i/a 

=  ^  |df  lH(f)l^  ^  exp(-i2irfAn)  g(nA)  A  .  (E-2) 

'/a  n 

upon  use  of  (68)  and  (B-6),  along  with  the  observation  that  b(n)  in  (62)  is 
even  in  n;  the  integrals  over  frequency  f  are  over  any  interval  of  length 
1/A.  But  the  summation  over  n  in  (E-2)  is  expressible  as 

Jdt  exp(-i2irfT)  g(T:)  ^  A  6(r-nA)  = 

n 

=  G(f )  ®  2  6(f  -  ^  G(f  -  ^^1)=  G(f )  .  (E-3) 

n  n 

which  is  the  aliased  version  of  G(f)  in  (E-1).  Here  we  used  the  fact  that  the 
Fourier  transform  of  a  product  is  the  convolution  of  the  corresponding  Fourier 
transforms.  Employment  of  (E-3)  in  (E-2)  yields 
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2  s'"'''  ■  bw 


df  |H(f)l  ^  ^(f) 

’U  * 

’/^  - 


(E-4) 


in  terms  of  the  (periodic)  power  transfer  function  of  the  filter  H  and  the 
aliased  spectrum  G. 


Now  identify  general  function  g  as 

g(r)  =  , 

in  which  case 

G(f)  =  P^'^^f)  =  exp(-i2irfr)  p^r)  , 

according  to  (47).  Then  (E-4)  becomes 


(E-5) 


(E-6) 


^  B(n)  p^na)  =  ^  J  df  lH(f)|  ^  T^^'^^f)  ,  (E-7) 

n  '/a 

(k) 

in  terms  of  the  aliased  spectral  versions  of  P  (f). 

If  the  spectral  width,  l/T,  of  lH(f)l  ^  is  narrow  relative  to  the 
bandwidth  W  of  the  spectral  functions  in  (E-7),  (see  figures  5  and  6),  then 
an  approximation  to  (E-7)  is  afforded  by  observing  that 

.sria 

2  B(n)  P^na)  =  ^  j'  df  )H(f)l  ^  P^'^^f)  = 
n  -.s/a 

.Sk 

=  ^^^*^^0)  J  df  lH(f)|  ^  i  P^'^^0)  ,  (E-8) 

-.S/4 

and  using  (B-4). 
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According  to  (E-3)  and  (E-6),  we  have,  for  k  =  1, 

P<')(f)  =  <E- 

n 

which  is  the  aliased  version  of  the  input  spectrum  to  the  detection  system  of 
interest. 


For  k  =  2,  observe  first  from  (E-6)  that 

p(2)(f)  ^  ®  p(''^(f)  =  J  du  P^^\u)  P^^^f-u)  .  (E-10) 

Then  the  aliased  version  is 

p(2)(f)  ,  2  p‘2>(f  -  '■‘'’(f  ■-  S  -  - 

n  n 

=  j"  du  P^^^u)  P^^^f-u)  =  r  du  P^^^(u)  P^^^(f-u)  ,  (E-11) 

'  -9f(  1  1 

by  use  of  (E-9)  and  the  periodicity  of  P'  .  This  last  relation  states  that 
•  2 ) 

the  aliased  spectrum  P'  '(f)  can  be  found  by  convolving  the  aliased  input 

-d ) 

spectrum,  P  (f),  over  one  period. 


In  particular,  a  special  case  of  (E-11)  is 


= 


using  the  even  character  of  the  input  spectrum. 


(E-12) 
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APPENDIX  F.  VARIOUS  INPUT  SPECTRA 


The  input  spectra  considered  in  this  appendix  are  characterized  by  shape 
for  f  in  the  range  and  are  symmetric  about  the  origin  f  =  0.  We 

still  define 


(F-1) 


However,  the  effective  bandwidth  of  the  positive  frequency  components  of  this 
spectral  shape  is 


W 


e 


df  f 


df  f 


2n 


2n  +  1  1^2  M  } 

,  ,.2  £2n+l  ,2n+l 

(n+1)  f-  -  f^ 


and  must  be  accounted  for,  in  the  evaluation  of  the  standard  output 
signal-to-noise  ratio  in  (51). 


(F-2) 


The  results  for  the  normalized  correlation  and  the  effective  bandwidth 
are  listed  below  for  various  values  of  n,  where  we  use  the  abbreviations 

<*1  =  2irf.|r  ,  ^2  =  2-af^T.  (F-3) 
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n  =  -2:  -6  dB/octave 


COS(o^)  COS(a2) 


12  Q  +  1 


n  =  -1 :  -3  dB/octave 


CKa^)  -  Ci(a^) 


■  in(a2/aT)  ’ 


P(T) 


sin(o2)  -  sin(a^ ) 
“2  ■  “1 


M 


e 


=  W  . 


n  =  1 :  +3  dB/octave 


a_  sin(a-) 
p(T)  =  2  ^ 


sin(a^)  +  C0S(a2)  -  COS(a^) 


2  2 
“2  ■  “1 


W 


e 


=  W 


12  0 

2 

12 


2 


+  1 


(F-4) 


(F-5) 


(F-6) 


(F-7) 
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n  =  2:  +6  dB/octave 


pCc) 


2)  sin(x)  +  2  X  cos(x) 


] 


144  +  24  0^  +  1 

144  +  72  +  1  .8 


(F-8) 


The  way  in  which  effective  bandwidth  enters  the  standard  output 
signal-to-noise  ratio  is  via  (51);  namely 


since 


[[df 

Jdf 


(F-9) 


(F-10) 


If  we  equate  (F-9)  to  (158),  and  solve  for  the  required  input  signal-to-noise 
ratio,  we  obtain 


S/N 

(S/N)^ 


(F-11) 


where  ^  is  given  by  (159);  the  last  equality  in  (F-11) 
-6  dB/octave  spectrum  and  has  utilized  (F-4). 


applies  only  to  the 
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APPENDIX  G.  EVALUATION  OF  (166)  FOR  FLAT  BANDPASS  SPECTRUM 

We  are  interested  here  in  calculating  the  quantity 

oo 

S(k)  =  (nA)  for  k  >  1  ,  (G-1 ) 

n=l 

where  the  relevant  normalized  correlation  is 

p(T)  =  cos(2irf(;T)  sinc(WT)  .  (G-2) 

Sampling  increment  A  is  arbitrary  in  relation  to  center  frequency  f^  and 
bandwidth  W. 

Frequency  Domain  Representation 

k  f  kl 

Since  p  (t)  and  P'  ■^(f)  are  related  by  Fourier  transform  (47), 

P^'^^f)  =  j  dT  exp(-i2^fT)  p^T)  .  (G-3) 

their  samples  are  related  through  Poisson's  formula 

1  /(n.)  =  1  i  ; 

n=--o  n=-<>4 

see  [5,  page  36,  (36)],  for  example.  This  relation  is  extremely  useful  for  a 
bandlimited  input  spectrum,  since  its  k-fold  convolution  P^*^^(f)  will  have 
limited  extent  in  f,  leading  to  finite  summations  on  the  right-hand  side  of 
(G-4). 
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A  more  appropriate  version  of  (G-4)  for  our  use,  where  p{z)  is  real  and 
even,  is  furnished  by  the  following: 


CO 


S(k)  =  ^  P^nA)  =  1 
n=l 


CO 


p^nA)  -  1 


_n=-«> 


oo 


1_ 

2A 


n=-<so 


n=l 


Evaluation  of  Auxiliary  Functions 


Before  directly  evaluating 


P^'^^f)  =  [  exp(-i2irfT)  [cos(2irf^r)  sinc(Wr^'^  , 


(G-6) 


we  first  consider  the  auxiliary  functions 


G^(u)  =  1  dx  exp(-i2.ux)  sinc^x)  for  k  >  1 


(G-7) 


It  is  readily  shown  that 


G^(u)  =  rect{u)  =^1/2  for  juj  =  1/2  >. 


for  )ul  <  1/2^ 
for  juj  =  1/2 
for  lu]  >  1/2 


(G-8) 


By  repeated  convolution,  there  then  follows  (see  [6,  pages  11,  12,  33,  34]) 


G2(u)  = 


.f 


-  |u|  for  luj  <  1 


0  for  luj  >  1 


I  -  ^ 


for  ]u|  <  2 


G  (u)  =  <  i(3  -  2lul  for  j  <  |u|  <  f 


for  juj  >  2  j 


(G-9) 


(G-10) 


106 


TR  7633 


1(4  -  6u^  +  3  lul^)  for  juj  <  1 


G.(u)  =<1(2  -\n\)^ 


for  1  <  lu)  <  2 

0  for  |u|  >  2 

Generally,  G|^(u)  =  0  for  ju)  >  k/2. 

Evaluation  of  p(*^)(f) 

By  expanding  the  k-th  power  in  (G-6)  according  to 


(G-11) 


cos'^(y)  =  Ij^  jjxp(iy)  +  exp(-iy)]*^  , 


(G-12) 


and  using  (G-7)-(G-11 ) ,  there  readily  follows 

1  p(i),n  _  ^  L  .  G 

A  ^  2WA  ^1\  W  j~  2WA  L^H  W  y  ^1  W  Jj’ 


1  p(2)(fs  =  _L 

A  ^  4WA 


1  p(3)(f)  =  ^ 

A  ^  ’  8WA 


^  3, 


i  p(4)(f)  ^  J_ 

A  ^  ’  16WA 


(f  +  4f 


4V  W 


+  4G, 


ff  ±  2f^ 
~W 


+  6G 


,(S) 


CG-13) 


(G-14) 


(G-15) 


(G-16) 


where  the  +  shorthand  notation  is  explained  in  (G-13) 


Equations  (G-13)-(G-16)  coupled  with  (G-8)-(G-ll)  enable  us  to  evaluate 
the  desired  quantity  in  (G-5)  for  k  =  1,2, 3, 4,  in  a  very  efficient  manner, 
since  all  the  sums  in  (G-5)  terminate  after  a  finite  number  of  terms.  On  the 
other  hand,  for  k  >  5,  we  resort  to  direct  numerical  evaluation  of  (G-1)  and 

5 

(G-2),  since  the  summands  decay  at  least  as  fast  as  1/n  and  can  be 
terminated  with  a  desi.red  level  of  error. 
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Error  Analysis 


Substitution  of  (G-2)  into  (G-1)  and  use  of  A  =  T/M  yields 


S(k)  =  ^  rcos^2irQ^  n^  sinc^^  n^ 

n=l  ^ 


(G-17) 


If  this  sum  is  conducted  through  N  terms,  the  error  E  is  upper  bounded 
according  to 


oO 


CO 


<-  1 


1 


1 


,  ,k  -  k  ) 

«  M.T  (o^n)  a  ^  X 

n=N+l  '  '  jij 


dx 

k 


1 


iL.k-l 

a  (k-1)  N 


(G-18) 


where  a  =  irTW/M.  Solving  for  the  N  required,  we  obtain 


N  >  N.  =  “  exp 


-inf(k-l)  aE] 


k-1 


(G-19) 


which  depends  on  k,  as  well  as  a  and  specified  error  E.  This  limit  was  used 
in  (G-17)  for  k  >  5,  with  E  =  l.E-10,  to  get  the  results  in  figures  33  and  34. 
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APPENDIX  H.  PROGRAM  LISTINGS 


In  this  appendix,  we  collect  all  the  programs  that  were  used  to  generate 
the  numerical  results  and  figures  in  the  main  body  of  the  report.  The 
auxiliary  functions  Si(x)  and  r(x)  common  to  these  programs  are  listed  in 
table  H-6. 


In  order  to  convert  table  H-1  to  the  Gaussian  spectrum  instead  of  the 
flat  bandpass  spectrum,  the  following  changes  and  additions  are  required: 


1  0  0 

Te=. 

1  1  1 

F=Te*Hs+Ns 

112 

IF  F>100.  THEN  210 

120 

Rho  =  COS<T.:^Ni:>^EXP  (-F) 

!  GRUSSIflN  SPECTRUM 

In  order  to  convert  table  H-2  to  the  Gaussian  spectrum,  the  following 
changes  and  additions  are  required: 


1  0  0 

Te  =  .  = 

li^^P 

I*Tui^Ti..i 

1  1  1 

F  =  Te7 

★  Ns 

1  12 

IF  f; 

>10 

0.  THEN 

2  0  0 

1  20 

R  h  o  ~  1. 

I:  III  S: 

(Tc^Ns) 

★  EXP 

-F> 

!  GAUSS  I  AH  SPE 

CTPUM 

irder 

to  convert 

table  H 

-1  to 

the  RC 

filter  instead  of 

the  box 

filter,  use: 

101  Te  =  2..ri 

110  FOP  Hs=l  TO  M*15 

140  PI:  =E:-':P -Te*Ms  >  !  PC:  FILTER 
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In  order  to  convert  table  H-2  to  the  RC  filter,  use: 

101  Ti==2. 

110  FOR  N£=1  to  N*15 

130  Pk=EXP';-T6*Hs>  •!  RC  FILTER 

Table  H-5  lists  the  program  used  for  investigating  the  long  averaging 
time  assumption;  however,  it  is  not  recommended  for  use,  since  the  earlier 
programs  are  capable  of  exact  evaluation  of  required  signal-to-noise  ratios. 

The  word  DOUBLE  denotes  INTEGER  variables  in  Hewlett-Packard  BASIC  on  the 
9000  computer. 
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Table  H-1 .  Full-Wave  Rectifier,  Flat  Bandpass  Spectrum,  Box  Car 


1 0  T  w  =  5  0 . 

20  Q=5 . /S . 

30  M=1000 

40  K=1000 

50  DOUBLE  M,K,Hs,Ks 
60  REDIM  S( 1:K) 

70  DIN  S<;i:1000> 

30  MRT  S='-'0.> 

30  T.:=2.  *PI*Q*Tw.-M 

100  Ts  =  PI*TW''’M 

110  FOR  Ns=l  TO  M-1 

1 2  0  R  h  o  =  C  0  S  (.  T  c  *  H  •=■  >  *  S  I  H  (  T  2-  *  N 

1 3  0  F;  h  o  2  =  R  h  o  *  R  h  o 

140  Pk  =  l.-H=./M 

150  F 0 R  K £  =  2  TO  K  STEP  2 


TW  FULL- WAVE  RECTI  FI 

fc/-W  >=  .5 

M ,  n  u  m  b  e  r  o  f  s  am  p  1  e  s 
N  umber  o  f  t  e  r  rn  s  i  n  St-i  m 


s  >  C  T s * H s  >  !  FLAT  BA H D P A S S  S P 

!  BOX.  CAR  FILTER 


160  Pk=Pk*Rho2 

170  S'::Ks;;'=sa(s>+Pk 

1  30  I  F  ABS  Pk  <  1  .  E-  1  2  THEN  200 

130  NEXT  Ks 

200  NEXT  Ns 

210  PR  I  NT  "  TW  =  "  ;  Tw ;  "  0.  = "  ;  Q ;  "  M  = "  ;  N ;  "  k 

220  INPUT  "NU  I  NIJ 


230 

G  5  =  F  G  -Sirn  rii  =l  k 

. 5*V+. 5) 

240 

U=1 . 

250 

3  u  rri  =  0  ■ 

2  6  0 

FOR  Ks=2  TO 

K  STEP  2 

270 

T=V+2. -Ks 

280 

U  =  IJ^T*T.--(Ks 

*  ( K  s  -  1  >  > 

2 '?  0 

Surn  =  Sun'i  +  U*S 

<  K  s.  > 

3  0  0 

NEXT  Ks 

310 

D  M  =  S  Q  R  •-!  F*  I  . '  ^ 

F  ^■■1  ij  ctrn  ru  -zL  V  +  .  5  )  •••*  U  5  Ij  5  .* 

3  2  0 

S  u  m  =  D  Kf  -  1  .  +  2 

.  ★Sufi'i 

330 

F=2. ★SQRCTw 

★  s « j  rii  •■■■  n  >  '•/ 

340 

PRINT  V,10. 

*  L  G  T  F  ,  2  .  LI  ^  S  '■!  K  )  ^  S-  ‘.4  rri 

350 

GOTO  220 

360 

END 

Fi Iter 

ER 

;CTRUM 
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Table  H-2.  Half-Wave  Rectifier, 

Flat  Bandpass  Spectrum,  Box  Car  Filter 

10 

Tw=50.  1 

TN  HRLF-NRVE  RECTIFIER 

20 

Q  =  5..--6.  1 

fc-'W  >=  .5 

30 

M=1000  ! 

^  n  u  m  b  €’  r  o  f  s  a  m  p  1 1-  s- 

40 

K=1000  ! 

N  u  m  b  e  r  'i'  f  t  e  r  m  •=■  in  Su  m 

50 

DOUBLE  M,K,H2.,Ks 

60 

REDIM  S'::i:K),IJ'::0:K) 

70 

D I N  S  <  1:1 0  0  0  > , U  <  0 ;  1 0  0  0  > 

8  0 

MRT  S=<0,> 

90 

Tc  =2 . I *Q*Tw/M 

1  00 

Ts  =  PI*Tw--M 

110 

FOR  Hs=l  TO  M-1 

120 

Rho  =  COS<Tc-i^H3>+SIN<Ts^Ns::' 

(Ts^Ns)  !  FLRT  BRNDPR33  3PECTRUM 

130 

Pk  =  l  .  -Ns/-M 

I  BOX  CRR  FILTER 

140 

FOR  Ks=l  TO  K 

150 

Pk=Pk*Rho 

160 

S<Ks>=S<Ks>+Pk 

1  70 

IF  RBS<Pk >< 1 . E-12  THEN  190 

130 

NEXT  Ks. 

190 

NEXT  Ns 

2  0  0 

PRINT  "TW  ="  ;  Ti..i;  “  Q  =";Q 

: ;  “  M  =  "  ;  M  ;  "  K  =  "  ;  K  ;  "  3  <  K  >  =  "  ;  3 

210 

INPUT  "NU 

!  NU 

220 

G5  =  FNGarrim.=i<  -  5*V+.  5> 

230 

LI  (  0  )  =  1  . 

240 

U  =  F  N  G  am  m  a  <  .  5  *  V  +  1  .  /■  G  5 

250 

U<  1  =2.  *U*U 

260 

S  u  m  =  U  <  1  >  *  S  <  1  > 

270 

FOR  Ks=2  TO  K 

230 

T=V+2. -Ks 

2  9  0 

U  <  Ks  =U  =  U  <  Ks-2  )  ★T-i^T.-'  <  Ks^  <  K 

:s-i  >  > 

300 

S  u  m  =  3 1  j  m  +  U  *  S  (  K  s  ') 

3  1  0 

NEXT  Ks 

3  2  0 

D  ••.«»  =  3  Q  R  (  F*  I  >  F  \\  G  am  m  a  < +  .  5  >  < 

5  ^  G  5 

330 

3  u  m  =  2 . +  D  M - 1 . +  2 . ^3  u m 

340 

F  =  2  .  *  3  Q  R  (  T  w  +  3  m  /  M  •  V 

350 

PRINT  V  ,  1  0  .  t^LGT  <  F  2  .  ^U^3  < 

K  )  •  •  :E:  i.j  rii 

3  6  0 

GOTO  210 

3  7  0 

END 

112 


TR  7633 


Table  H-3.  Full-Wave  Rectifier,  -6  dB/Octave  Spectrum,  Box  Car  Filter 


10 

Tw=50.  !  TW 

FULL-WAVE  RECTIFIER 

20 

Q  =  5..-"6-  !  fc-’-W 

30 

M=100  !  M,  number 

of  s 

amp  1 es 

40 

K  =  1 0  0  0  !  H  u  m  ti  €■  r  o  f 

t  e  r  m 

s  i n  Sum 

50 

DOUBLE  r‘1,K,Nsj,Ks 

60 

REDIM  S'::i:K> 

70 

DIM  1 :  1000;:' 

CO 

MAT  S='::0.::' 

30 

Cl=2.  *PI*Tw*<Q-.  5>.-'M 

100 

C2=2. *PI*Tw*<Q+- 5>/M 

110 

Sci=  1 2 . * 

120 

5  c^  =  5  U  R  S  3 .  )  S  c^  +  1  •  ) ) 

130 

FOR  Ns.=  l  TO  M-1 

140 

Rho  =  FNRho'::Hs,  Cl  ,  02 >  ! 

—  6 

dB. -'OCTAVE  SPECTRUM 

150 

R  h  o  2  =  R  ^'1  o  ^  R  h  o 

160 

Pk  =  l  .  -Hs.-'M  ! 

BO 

X  CAR  FILTER 

170 

FOR  Ks=2  TO  K  STEP  2 

. 

1S0 

Pk=Pk*Rho2 

190 

S<Ks>=S<Ks::'+Pk 

2  0  0 

IF  ABS<Pk  )<  1  .  E-r2  THEN  220 

210 

NEXT  Ks 

220 

NEXT  Ns 

230 

PR  I  NT  "  TW  =  '*  ;  Tw  ;  "  Q  =  "  ;  Q  ;  "  M  =  "  ; 

M ;  " 

K  =  "  ;  K  ;  "  S  K  >  =  “ 

240 

INPUT  "NU  =%'•/  ! 

NU 

250 

G5  =  FNGamma‘::  .  5^V+.  5> 

260 

U=1  . 

270 

3  u  m  =  0  • 

230 

FOR  Ks=2  TO  K  STEP  2 

290 

T=V+2. “Ks 

300 

ijzilj^T^T.-  ^Ks^kKs-l  )  > 

3  1  0 

S  u  m  =  S  u  m  +  U  *  S  <  K  ■=■  ) 

3  2  0 

NEXT  Ks 

3  3  0 

D  M  =  S  Q  R  <  F’  I  >  ★  F  N  G  arii  rri  a V  +  .  5  >  ■  ( i3  5  *  i3  5  ) 

340 

S  u  m  =  D 1  .  +  2  .  ^  S  u  m 

350 

F  =  2  .  ^  S  C!  R  (  T  w  ^  S  u  rn  ■  [*1  > .  ■  V  *  S 

3  6  0 

P  R  I  N  T  V  ,  1  0  .  -fr  L  T  (  F  )  ,  2  .  ^  l-l  ^  'S  (  K  • '  '3  u  rri 

370 

G 0 TO  24 0 

3  3  0 

END 

390 

1 

4  0  0 

DEF  FNRho'.  DOUBLE  Ns,  REAL  C1,C2::' 

410 

Al=Cl*Ns 

420 

R2=C2*Ns 

4  3  0 

R  h  o  =  F  S  i  A  1  >  -  F  N  3:  i  (32  •'  +  C-  Ci  '3  (  A  1  >  A  1  -  i3  0 :5 

:  ( A2 

.--A  2 

440 

R  h  o  =  A  1  *  A  2  F!  kj  C‘  A  2  “  A  1 

4  50 

RETURN  Rho 

460 

FNEND 
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Table  H-4.  Half-Wave  Rectifier,  -6  d&/Octave  Spectrum,  Box  Car  Filter 


10 

T  w  =  50 . 

j 

TW 

HflLF-NFiVE  RECTIFIER 

20 

Q  =  5 .  ••  •*’  6 . 

1 

fc.-W 

30 

M=100 

j 

t'1  J  n  u  rn  b  s-  r 

o  f  s  ari'i  p  1  €'  S' 

40 

K=1000 

! 

h-l  u  rn  b  6-  r  o  f ‘ 

t  e  r  rn  s  i  n  :=:  u  rn 

50 

DOUBLE  M,K,Hs.,Ks 

60 

REDIM  S<  1  s  K>  ,  IJ^O:  K> 

70 

D I M  S  1 :  1 0  0  0  > ,  U  0 :  1 0  0  0 

> 

CO 

MAT  S=<0.::' 

30 

C  1  =  2  .  *  P  I  *  T  I...I  *  C  Q  -  .  5  >  ■  M 

100 

C2  =  2.  5>/M 

1  10 

Sq=12. *Q*Q 

120 

S  =  S  Q  R  (  (  S  c^  —  3  .  )  .•  *’  (  S  +  1  .  ) 

‘]i 

130 

FOR  Ns=l  TO  N-1 

140 

Rho=FHRho<Ns, Cl , 02 > 

!  -6  dB.'-OCTRVE  SPECTRUM 

150 

Pk  =  l .  “Hs.-  M 

!  BOX  CRR  FILTER 

160 

FOR  Ks=l  TO  K 

170 

Pk=PMRho 

1S0 

SCKs>=Sa<s>+Pk 

190 

IF  RBS<Pk >< 1 . E-12  THEN 

2 

10 

200 

NEXT  Ks 

210 

NEXT  Ns 

220 

PRINT  "TW  =";Tw;‘'  Q 

=  '« 

;  Q  ;  ■*  M  =  " 

;  M  ;  K  =  "  ;  K  ;  “  S  (  K  >  =  '•  ;  S  <  k 

230 

INPUT  "NU  =",V 

!  NU 

240  • 

G5  =  FNG.arrirri-=t(  .  5*V+.  5> 

250 

U  (  0  >  =  1  , 

260 

U  =  FNG.=Lrrirri.=L(  .  5*V+  1  .  >  -"GZ 

270 

U'::  1  >=2.  *U^U 

2  S  0 

S 14  hi  =  U  <  1  >  *  S  <  1  > 

2  9  0 

FOR  Ks=2  TO  K 

3  0  0 

T=V+2. -Ks 

3  1  0 

U  <  Ks  >  =U  =  U  Ks-2  ^T^T--^  (  K 

S^ 

Ks-1  > 

, 

3  2  0 

S  u  h  i  =  S  u  rn  +  U  *  S  <  K  s  > 

330 

NEXT  Ks 

340 

D  M  =  S  Q  R  <  F'  I  >  ^  F  h-l  G  arn  rri  a  <  +  . 

Cj 

G  5  ^  i3  5 

3  5  0 

S  u  m  =  2  .  ★  D  -  1 .  +  2  ,  S  L4  ro 

3  6  0 

F  =  2  .  Q  R  T  i.j  ^  1,4  ri'i  ■  '■/  ^  S 

q 

370 

PRINT  V,  10. +LGT-::f:>  ,  2.  ^ 

u* 

E:  C  K  !:*  /  S  i.i  rri 

3  3  0 

GOTO  230 

3  9  0 

END 
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Table  H-5.  Half-Wave  Rectifier,  Flat  Bandpass  Spectrum, 


10 

Long  Averaging  Time  Assumption 

Tw=50.  !  TW 

HflLF-WflVE  RECTIFIER 

20 

Q=. 501 

!  fc/W  >=  .5 

LONG  RVERflGING  TIME 

30 

M= 1000 

!  M  j  number  of  s-arn 

pies  ASSUMPTION 

40 

K= 1 000 

!  H  u  m  b  e  r  o  f ‘  t  e  r  m  s 

i  n  Sum 

50 

60 

RED  I M  SO:  K  > ,  U  (  0 :  K  > ,  H  <  5 :  K  ;■ 

D I M  S  <  1:1 0  0  0  > , U  <  0 :  10  0  0  > , R  h  o  <1:1 0  0  0  0  > , R  k  <  1  : 

1  0  0  0  0  !:> 

70 

DOUBLE  N*::5:  10000::- ,  M, 

K, Ns, Ke, Nmax 

80  Q2=Q+Q 

90  Q3=Q2+Q 

100  Q4=Q3+Q 

110  S1=FNG1  Olf::' 

120  S  2  =  F  N  G  2  <  Q  2  >  +  F  N  G  2  <  0 . > 

1  3  0  S  3  =  F  ^■'1 G  3  (  Q  3  )  +  3  .  *  F  f  ^  G  3  '•!  Q  ) 

140  S  4  =  F  N  U  4  L-!  4  !•'  +  4  ,  ^  F  h-l  G  4  k  i3  2  )  3  .  N  i  j  4  (  0  .  ) 

1 50  Owd  =  M--'Tu 

160  F  o  w  =  0 • 

170  F  o  w  =  F  o  w  +  0  w  d 

130  IF  Fow-Q4>=2.  THEN  250 

1  9  0  'z>  1  =  1  +  F  h-l  G  1  F  o  ui  ~  Q  )  +  F  h-l  G  1  (  F  o  w  +  il!  } 

200  S  2  =  S  2  +  F  N  G  2  <  F  o  w  -  Q  2  )  +  2  .  ^  F  N  G  2  F  o  w  >  +  F  N  G  2  F  o  m  +  Q  2  ) 

2  1  0  ‘6  3  =  3  3  +  F  I j  3  F  c>  w  “  Q.  3  )  +  - 5  •  N  G  '3  •'!  F  o  i.*j  - 1“!  )  +  3  .  *  F  h-l  i Ij  3  (  F  '3  i.'j  +  0.  )  +  F  h-l  G  3  (  F  C'  «.>.»  +  Q  3  ) 

220  S4  =  S4  +  FMG4(Fqw-Q4::‘+4.  ^FNG4  <  FciUi-Q2  > +6  ,  ^FNG4  Foi.j  > 

2  3  0  3  4  =  34  +  4.  +  F  h-l  G  4  '•!  F  o  1.0  +  Q  2  !•’  +  F  t  ^  G  4  '■!  F  o  j.'J  +  Q  4  ) 

240  GOTO  170 

250  S  < 1 >  = . 5  +  0  w  d  *  S 1 - . 5 
260  3  <  2  >  =  .  25^0i...id*S2~  .  5 

2  7  0  S  (  3  >  =  ,  1  2  5  *  Cl  I  j  id  +  S  3-  .  5 

2  3  0  S  (  4  )  =  .  0625  +  0  ui  d  +  C;  4  -  .  5 

290  T  c  =  2  .  +  P  I  +  Q  +  T  IJ  N 

300  Ts  =  PI+Tij-'M  ‘ 

3  1 0  E  r  r  o  r-  =  1  .  E  -  1  0 

320  Cl=l./Ts 

3  3  0  C  2  =  -  L  0  G  K  T  s-  +  E  r  r  o  r  > 


340 

FOR  Ks=5  TO  K 

350 

T=Ks«l 

3  6  0 

N  (  K  s  >  =  I  N  T  (  C  1  +  E  X  P  <  C  2  -  L  0  G  T  >  '  T 

'  ;:-  + 1 

370 

NEXT  Ks 

330 

FOR  Ks=K-l  TO  5  STEP  -1 

3  9  0 

IF  N  K  s  <  N  •:  K  s  +  1  THEN  N  (  K  s  =  N  h 

:;s+l  > 

400 

NEXT  Ks 

4  1  0 

N  m  h.x  =  M  A  X  N  +  .-  > 

420 

R  E  Ii  I  M  F!  Li  o  (  1  :  N  m  a;-:;  ,  F!  k  1  :  N  m  a> 

430 

FOR  Ns=l  TO  Nrria::< 

440 

F:  =  Rho  (  Ns  =C0S  Tc  +Ns  >  +S  I  N  (  Ts  +  Ns  :: 

■••■■(Ts  +  Ns 

4  5  0 

Rk  Ns:;'=F:+R+R  +  F: 

4  6  0 

NEXT  Ns 

4  70 

FOR  Ks=5  TO  K 

430 

S  =  0. 

4  9  0 

FOR  Ns=l  TO  N'::Ks::' 

J  0  0 

R  k  =  R  k  N  s  ‘  =  P  k  (  N  s  '  ^  F:  L.  <3  ■  ‘ '  s  ;• 
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Table  H-5 


ont'd).  Half-Wave  Rectifier,  Flat  Bandpass  Spectrum, 
Long  Averaging  Time  Assumption 


510 

S=S+Rk 

520 

NEXT  Ns 

5:30 

S<Ks>=S 

540 

NEXT  Ks 

550 

PRINT  “T 

560 

INPUT  "N 

570 

G  5  =  F  N  G  arr 

580 

U  <  0  )  =  1  . 

5  9  0 

U  =  FNGariUi 

6  0  0 

LI  C  1  >  =2  .  ^ 

610 

S  u  rrt  =  U  ( 1 

620 

FOR  Ks=. 

630 

T  =  V  +  2.  -I 

640 

U‘;Ks>=U 

650 

S  u  ffi  =  S  u  rn 

6  6  0 

NEXT  Ks 

670 

D  V  =  S  Q  R 

680 

Sum=2 . * 

6  9  0 

F=2. *SQ 

-700 

PRINT  V 

710 

GOTO  5t 

720 

END 

730 

1 

740 

DEF  FNG 

750 

fl=fiES<F 

760 

IF  fl<.^ 

770 

IF  fl=." 

780 

.RETURN 

790 

FNEND 

300 

1 

3  1  0 

DEF  FN 

3  2  0 

fl  =  flBS‘:: 

330 

IF  fi<l 

340 

RETURN 

850 

FNEND 

3  6  0 

I 

370 

DEF  FN 

380 

fl=fiES( 

390 

IF  Fl<. 

9  vj  0 

IF  fl<l 

910 

RETURt 

9  2  0 

FNEND 

930 

1 

9  4  0 

DEF  Ft 

9  5  0 

fl  =  flE3 

960 

fl2  =  H*t 

970 

IF  H< 

9  3  0 

IF  Fl< 

9  9  0 

RETUR 

1  0  0  0 

FNEtID 

=";TU!"  Q 

—  II  ',,1 
“  S  ■ 

, 5*V+. 5> 

.  5*V+1  ■ 


M 


K  = " ;  K ; 
NU 


S'::  1 ) 

TO  K 

i*S'::Ks> 


.  >*FNGarnma'::V+.  5  > G5*G5  > 
- 1 .  +  2  -  u  rn 
T  1.0  *  S  1.4  rft  • 

I  0  .  *  L  G  T  <  F  >  ,  2  .  *  U  *  S  K  )  S  u  ni 


<  F  > 

THEN  RETURN  1. 
THEN  RETURN  .5 


AF> 


THEN  RETURN  1 . -H 


3-::F> 

I 

THEN  RETURN  . 75-fl*fl 
S  then  return  .  125*':3. -fl-fl  j 


4(F> 


THEN  RETURN  (  4  . -6  .  *Fl2  +  3  .  *fi2^H  ^  _  b 
THEN  RETURN  ■.  2  . -H  j  *  (  2  . -H  '  ‘72  .  H  ,■ 

0  . 


S<K) 


S  K 
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Table  H-6.  Auxiliary  Functions  Si(x)  and  r(x) 


10 

DEF  FNS  i  < X  >  !  i  ( jo  y  i  a  5 . 2 .  S ,  1 4 , 3S ,  39 

Y  =  X*K 

20 

30 

IF  Y>1.  THEN  70 

40 

i  •  y  •-»  4  4  6  1*  1  ci  y  1  8  E ""  5  “  T'  ,  0  6  1  9  2  4  3  5  8  2  2  E  “  7 

50 

Si  1  .  .  0555555555556-V*.::  1 . 66666666667E-3- Y 

^  l"l 

i  )  >  > 

60 

RETURN  Si 

70 

7 1  —  •-*  o  •  1  y  c-  4  ?  5  +  V  3  35. 6  r'  7  3  3  y  +  Y  1^!  6  5  ■  1 8  7  0  3  3  +  y'  *  '•!  Si  8 . 0 

Cm  r* 

264  +  Y ) ) ) 

8  0 

T  2  =  1  5  7  .  1  0  5  4  2  3  +  Y  *  <  570. 2  3  6  2  8  0  +  Y  ^  Si  2  2 . 62491  1  +  Y  *  4  0 

02 

1433  +  Y>  > 

90 

F  =  T1.-'':;X*T2;:' 

1  0  0 

Tl=21 . 82 1 899  +  Y* < 352 . 0 1 8498  +  Y*  C302 . 757865  + Y* < 42 . 2 

42 

855+ Y) ) > 

110 

T2  =  449. 690326  +  Y+':;  11  14. 978885  + Y+ 482 . 485984  + Y+ 49 

■  196927  +  Y ) 

120 

G  =  T1..-.:;y*T2:> 

130 

S  i  =  1 . 57079632679  +  SGN  X  >  -F  +  COS  (  X  >  -G  +  9  I  N  y'  X 

140 

RETURN  Si 

150 

FNEND 

10 

DEF  F  N  G  .=Lni  rn  a  (  X  > 

!  HRRT,  page-  135,  #5243 

20 

DOUBLE  N,K 

3  0 

N=INT(X) 

40 

R  =  X~N 

50 

IF  N>0  OR  R>0.  THEN  80 

60 

PRINT  ■'  F  N  G  am  m  a  (  X  >  I  S  N  0  T 

DEFINED  FOR  X  =  ";X 

70 

STOP 

8  0 

^IF  R>0.  THEN  110 

90 

b  a  nun  a2=  1  . 

1  0  0 

GOTO  180 

1  10 

P  =  4  39 . 3 30 444 06 00 256 7 6 +  F:+  >: 

:  50 .  1 086937529709530 +  R  + 6 . 

120 

P  =  8762 . 7  1  0297852  1  4896  +  F:+  >: 

:2008. 527401  307279  12  +  R  +  P 

1  30 

P  =  42  353 . 6S95097440896  +  F;+  •: 

:  20886 . 86 1 7892698874  +  R  +  P  > 

140 

0  =  499. 02852662  1  439048-R+  •: 

.  1  i:i  9 . 4  9  Si  2  3i  4  1  5  7  01 2 13  0*  1  6  -  F!  +  (  2 

150 

Q  =  994  0  .  Si  07415  0  8277  2i  9  Si  -  F:  +  ■: 

.  1  5  2 13 . 6  01 7  2  7  Si  7  7  9  5  2  2  01 2  +  F!  +  0!  > 

1  60 

0  =  4  2353 . 6895097440900  +  F:  +  >: 

:  2  9  S:  0  .  Si  S  5  S:  S:  0  9  2  5664  9  9  -  R  + 1.-! 

170 

ij  am  rn  a2  =  P  0 

!  b  arri  rrj  a  2  +  P:  f  izi  r-  0 

1  80 

IF  N>2  THEN  220 

1  90 

IF  N<2  THEN  270 

2  0  0 

b  a  rn  r  1 1  a  =  b  a  rn  rn  a  2  ^ 

2  1  0 

Pi  E  T  U  P'  f  i  b  arn  rn  a 

2  2  tt 

b  arn  rn  a=  G  arn  rn  a 2 

2  3  0 

FOR  K=1  TO  N~2 

240 

b  aril  rn  a= b  arn  rn  a+  X  “  K  ) 

250 

NEXT  K 

260 

RETURN  Gamma 

270 

P:=l . 

2  8  0 

FOP:  K  =  0  TO  1-N 

2  9  0 

R  =  P:  ^  +  K  > 

S:  0  0 

NEXT  K 

>:  1  0 

b  arn  ru  a = b  arn  ni  a2  i*' 

3  2  0 

RETURN  G  anirna 

3  3  0 

FNEND 

44950724592 


. 0S1551 


245 


1 


R  > 
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